Quasi-linear regime and rare-event tails of decaying Burgers turbulence 
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We study the decaying Burgers dynamics in d dimensions for random Gaussian initial conditions. 
We focus on power-law initial energy spectra, such that the system shows a self-similar evolution. 
This is the case of interest for the "adhesion model' in cosmology and a standard framework for 
"decaying Burgers turbulence" . We briefly describe how the system can be studied through pertur- 
bative expansions at early time or large scale (quasi-linear regime). Next, we develop a saddle-point 
method, based on spherical instantons, that allows to obtain the asymptotic probability distribu- 
tions Vlrjr) and VlQr), of the density and velocity increment over spherical cells, reached in the 
quasi-linear regime. Finally, we show how this approach can be extended to take into account the 
formation of shocks and we derive the rare-event tails of these probability distributions, at any finite 
time and scale. This also gives the high-mass tail of the mass function of point-like singularities 
(shocks in the one dimensional case). 

PACS numbers: Valid PACS appear here 
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I. INTRODUCTION 

The Burgers equation [U, [3, [1], which describes the 
evolution of a compressible pressureless fluid, with a non- 
zero viscosity, was first introduced as a simphfied model 
of fluid turbulence, as it shares the same hydrodynami- 
cal (advective) nonlinearity and several conservation laws 
with the Navier-Stokes equation. It also displays strong 
intermittency, associated with anomalous scaling expo- 
nents for the velocity structure functions, but this arises 
from the formation of shocks (i.e. singular structures in 
the inviscid limit v 0^) where energy is dissipated, 
whereas the structures that appear in Navier-Stokes tur- 
bulence seem to be more varied and less singular (be- 
cause of pressure effects) Nevertheless, due to its 
greater simplicity - it can actually be explicitly integrated 
through the Hopf-Cole transformation - the Burgers 
dynamics retains much interest for hydrodynamical stud- 
ies, particularly as a useful benchmark for approxima- 
tion schemes On the other hand, the Burgers equa- 
tion also appears in many physical problems, such as the 
propagation of nonlinear acoustic waves in non-dispersive 
media Q, the study of disordered systems and pinned 
manifolds or the formation of large-scale structures 
in cosmology (Tol . [Tl| , see ^] for a recent review. In the 
cosmological context, where one considers the inviscid 
limit without external forcing, it is known as the "ad- 
hesion model" and it provides a good description of the 
large-scale filamentary structure of the cosmic web [l^ . 
Then, one is interested in the statistical properties of the 
dynamics, as described by the density and velocity fields, 
starting with a random Gaussian initial velocity [1, [l^ 
and a uniform density. These initial conditions are the 
signature of quantum fluctuations generated in the pri- 
mordial Universe and agree with the small Gaussian fluc- 
tuations observed on the cosmic microwave background. 
In the hydrodynamical context, this setup corresponds 
to "decaying Burgers turbulence" [isj . 

This problem has led to many studies, focusing on 



power-law initial energy spectra (fractional Brownian 
motion) in one dimension, E(){k) oc fc", especially for the 
two peculiar cases of white-noise initial velocity (n = 0) 
l|,l 2|.ll4l . ll5p or Brownian motion initial velocity {n = —2) 
311a [l3, [3 • Indeed, in these two cases the initial 
velocity field is built from a white-noise stochastic field 
(either directly or through one integration), which gives 
rise to Markovian processes and allows to derive many 
explicit analytical results. For more general n, it is not 
possible to obtain full explicit solutions, but several prop- 
erties of the dynamics are already known 8, 13]. In par- 
ticular, for —3 < n < 1, the system shows a self-similar 
evolution as shocks merge to form increasingly massive 
objects separated by a typical length, L{t) - the integral 
scale of turbulence - that grows as L{t) ~ i2/(n+3)^ while 
the shock mass function scales as ln[n(> m)] ^ —rn^^^ 
at large masses [H, [13, IH, HO] • In spite of these common 
scalings, the range —3 < n < 1 can be further split into 
two classes, as shocks are dense for — 3 < n < — 1 but 
isolated for -1 < n < 1 p^ . 

In this article, we consider the decaying Burgers dy- 
namics in d dimensions, for random Gaussian initial con- 
ditions and power-law initial energy spectra such that the 
system displays a self-similar evolution. This is in par- 
ticular the case of interest in the cosmological context, 
which shows a hierarchical evolution as increasingly large 
scales turn nonlinear as time goes on. Applying to the 
Burgers dynamics methods that have been used to study 
the coUisionless gravitational dynamics encountered in 
cosmology, we present a saddle-point approximation (in- 
stanton technique) that allows to derive some properties 
of the velocity and density fields in two regimes, i) the 
quasi-linear regime associated with early times or large 
scales, and ii) the rare-event tails of the velocity and den- 
sity distributions at any time or scale. 

This article is organized as follows. We first intro- 
duce in section [TT] the equations of motion and the initial 
conditions that define our system and we recall the geo- 
metrical interpretation of the Hopf-Cole solution of the 
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dynamics. We also define the overdensity, Ty^., and the 
velocity divergence (i.e. spherical velocity increment), 
0,., within spherical cells of radius r, that are the two 
quantities that we study in this paper. Then, we briefly 
describe in section Hill how the dynamics can be studied 
through perturbative expansions, that hold at early times 
or large scales, and we make the connection with the Zel- 
dovich dynamics that is equivalent from a perturbative 
point of view. Next, we present in section HVl a saddle- 
point approximation that allows to derive the asymptotic 
probability distributions V{rjr) and 7^(0^) reached in the 
quasi-linear limit (i.e. at early times or large scales). 
Then, we show in section fVl how to modify this approach 
to take into account shocks, and we derive the rare-event 
tails of these probability distributions, at any fixed time 
and scale. This also yields the high-mass tail of the mass 
function of point-like objects (shocks in the one dimen- 
sional case). Finally, we conclude in section IVTl 



II. BURGERS DYNAMICS 

A. Equations of motion and initial conditions 

We consider the rf-dimensional Burgers equation in the 
inviscid limit (with d > 1), 



dtU + (u.V)u = i/Au, 



0+, 



(1) 



for the velocity field u(x, i), and the evolution of the 
density field p(x, t) generated by this dynamics, starting 
from a uniform density po at the initial time t = 0. The 
latter obeys the usual continuity equation 



dtp + V.(pu) = and p(x, 0) = po. 



(2) 



Then, since there is no external forcing in Eqs.(IT])-(l2]), 
the stochasticity arises from the random initial veloc- 
ity Uo(x), which we take to be Gaussian and isotropic, 
whence (u) = by symmetry. Moreover, as is well-known 
3, if the initial velocity is potential, Uq = — V^qj it re- 
mains so forever, so that the velocity field is fully de- 
fined by its potential t), or by its divergence 6'(x, t), 
through 



u = -VV', 6' = -V.u = AV'. 
Normalizing Fourier transforms as 



dk e 



ik.x 



(3) 



(4) 



the initial divergence is taken as Gaussian, homoge- 
neous and isotropic, so that it is fully described by its 
power spectrum Pg^ (fc) with 

(^o) - 0, (^o(ki)^o(k2)) = (5i3(ki +k2)Peo(fci), (5) 

where we note Sd the Dirac distribution. In this article 
we focus on the power-law initial power spectra. 



Thus, the initial conditions obey the scaling laws 

A > : eo{x-'k) 'i?' A''-^ eoik), (7) 

0o(Ax) X-^ 0o(x), (8) 

where means that both sides have the same statis- 

tical properties. This means that there is no preferred 
scale in the system and the Burgers dynamics will gen- 
erate a self-similar evolution for —3 < n < 1, as seen 
in section IIIDI This is why we only consider the range 
—3 < n < 1 in this article. For the initial velocity and 
potential this yields for any A > 0, 

uo(Ax) A"^ uo(x), ^o(Ax) '= A^ "00 (x). (9) 

Since we have u(k, i) = i(k/fc^)(?(k, i), the initial energy 
spectrum is a power law, 

(uo(ki).Uo(k2)) = <5l5(ki + k2)Eo{ki), (10) 

with 

Eo{k)=k-^Pg,{k)ock"+^-''. (11) 

The initial velocity correlation at distance x reads as 



(uo(xi).Uo(x2)) = J dke^-^Eo{k) 



= (2^)^ 



' J^^iikx) 
dkk'^-^ ' \_ Eo{k) 
[kx) 2 -'- 



(12) 



where x = X2 — xi and Jd_^{kx) is the Bessel function 
of the first kind of order d/2— 1, whereas the initial one- 
point variance is 



(|uo|2) = / dkEo{k) - 



27r2 
IW) Jo 



dk k'^-^Eoik). 



Pg„ (fc) CX k 



n+3— 



with 



3 < n < 1. 



(6) 



(13) 

Thus, for — 1 < n < 1 the initial velocity correlation de- 
creases at large distance as the power law ([T2|) . in agree- 
ment with the scaling ([5]), while the one-point variance 
at a; = 0, Ea. (fT3| . diverges because of the contribution 
from high wavenumbers. Then, the initial velocity field 
is singular (e.g., a white noise for d — \ and n — Q) but 
this ultraviolet divergence is regularized as soon as t > 
by the infinitesimal viscosity For —3 < n < — 1 the 
integral (|13p shows an infrared divergence. In this case, 
the initial velocity field is no longer homogeneous and 
only has homogeneous increments (but the divergence 
is still homogeneous) [3] . Then, to build the initial veloc- 
ity from its divergence one must choose a reference point, 
such as the origin xq = 0, with Uo(xo) — 0, and define 
the initial velocity in real space as 

Uo(x) = y"dk(e''^''-e*''«)iio(k), for -3<n<-l. 

(14) 

Then, Equation (fT2l) no longer applies but the ini- 
tial second-order structure function, (|uo(x) — Uo(xo)P), 
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grows as x~^~^ . Note that because of the nonhnear ad- 
vective term in the Burgers equation ([1]) , the increments 
of the velocity field are no longer homogeneous for t > 0, 
which also means that the divergence 6'(x, t) is no longer 
homogeneous either. However, at large distance from 
the reference point (i.e. taking the limit |xo| oo or 
|x| oo), we can expect to recover an homogeneous sys- 
tem (in terms of velocity increments and matter distri- 
bution), see [2l[ for more detailed discussions. This can 
be shown explicitly for the case d = 1 and n = —2, where 
the initial velocity field is a Brownian motion [T7l.[l8j. On 
the other hand, we may add a low-A: cutoff A to the ini- 
tial power spectrum and restrict ourselves to finite times 
and scales where the influence of the infrared cutoff is 
expected to vanish for equal-time statistics. 



B. Density contrast and linear mode 

In order to follow the evolution of the matter distribu- 
tion we define the density contrast, (5(x, i), by 



<5(x,i) 



p(x, t) - po 



(15) 



Then, if we linearize the equations of motion (HJ-lHl) we 
obtain the solution 



0L(k,i) = 0o(k)e-'"= *, 5i(k,t) = 0o(k) 



1 



(16) 

where the subscript L stands for the "linear" mode. In 
the inviscid limit, — > 0+, this yields 

0i(k,<) -^o(k), 5L(k,t) =t0~o(k), (17) 

which could also be obtained by setting = in Eq.(IT]). 
Then, when we study the system at a finite time t > 0, 
we can as well define the initial conditions by the linear 
density field ^^(x, i), which is Gaussian, homogeneous 
and isotropic, with a power spectrum 

- 3 < ?i< 1 : Ps^ (fc, t) = f Pg, (k) cx (18) 

and an equal-time two-point correlation 



C5^(xi,X2) = (^L(xi,t)(5L(x2,t)) 

Jd -I (kx) 



(27r)2 / dkk" 
Jo 



(kx) 



Pi j/fc) cx (19) 



where x = |x2 — xi|. Note that for any n > —3 the 
initial density field is homogeneous, even though the ini- 
tial velocity only shows homogeneous increments when 
-3 < n< -1. 

Here we may add a few comments on the initial condi- 
tions that are relevant to the cosmological context. Let 
us first briefly recall how the Burgers equation (jT]) arises 
in this case. In the standard cold dark matter scenario 
[221, about 83% of the matter content of the Universe is 



in the form of a cold dark matter component, whereas or- 
dinary baryonic matter only forms the remaining 17% (in 
addition, there is a dark energy component, which is con- 
sistent with a cosmological constant in the Einstein equa- 
tions, which makes about 72% of the energy content of 
the Universe, while the previous two matter components 
only form the remaining 28%), see [1^. The cold dark 
matter has a negligible velocity dispersion (whence the 
label "cold" ) and it has only very weak non-gravitational 
interactions (whence the label "dark" , as it has only been 
"seen" through its gravitational effects so far). Then, it is 
well described as a pressure-less fluid coupled to its own 
gravity (here we focus on the late Universe, after the 
end of the radiation-dominated era, about 5 x 10'* years 
after the Big Bang, and on scales smaller than the Hub- 
ble scale, where the Newtonian approximation is valid). 
Therefore, the growth of matter density fluctuations is 
governed by the pressure-less Euler equation and the con- 
tinuity equation, coupled to the Poisson equation, in an 
expanding background [24| . Since the gravitational force 
derives from the scalar gravitational potential, it does 
not generate any vorticity, and any primordial vorticity 
is diluted by the expansion of the Universe (this only 
holds in the linear regime, as shell-crossings can generate 
vorticity in a non-perturbative fashion). Then, using a 
rescaling of time and velocity field, that brings out the de- 
viations from the mean Hubble flow (and also absorbs the 
effect of the uniform cosmological constant), and making 
the approximation that the velocity and gravitational po- 
tentials remain equal (this is exact in the linear regime 
and in one dimension, d — 1, before shell-crossing), one 
obtains the Zeldovich equation [25|. This corresponds 
to the Burgers equation (jT]) with i/ — 0. Then, one 
adds an infinitesimal viscosity, — > 0+ , to prevent shell- 
crossing [l3, [iH. This induces a sticking of particles 
within shocks, that is intended to mimic the trapping 
within gravitational potential wells [T^ . 

Next, in the cosmological context, the present mat- 
ter density fluctuations are assumed to arise from the 
growth of tiny quantum fluctuations generated during an 
inflationary stage in the early Universe. Moreover, these 
Gaussian initial fluctuations almost have a Harrison- 
Zeldovich power spectrum, that corresponds to n = 1 in 
Eq.p^ above (observations give n ~ 0.96 [1^). The case 
n = 1 is also called "scale-invariant" , as it gives a gravi- 
tational potential power spectrum of the form P^^ (fc) oc 
l^n-i-d _ f^-d^ gp ^^g^^ wavenumbers contribute with 
the same weight and the two-point correlation is formally 
scale- invariant, C^/,o(x) cx / dkk'^~^ PTf,g{k)W (kx) is inde- 
pendent of x (where W{kx) is some filtering function on 
scale x). Within the infiationary scenario, this property 
arises from the fact that the only relevant scale is the 
Hubble scale, that remains roughly constant during this 
stage (this can also be understood from the fact that 
during an exponential expansion there is no genuine ori- 
gin of time, i.e. the de Sitter spacetime is invariant un- 
der time translations, so that wavelengths generated at 
different times share the same properties). Then, since 
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these fluctuations have remained small until recent times 
they have evolved through linear theory until the matter- 
dominated era and the Newtonian regime. Therefore, 
they have remained Gaussian and different wavenum- 
bers have evolved independently (the linearized equa- 
tions of motion are diagonal in Fourier space) until a 
redshift z 10^. However, the primordial spectrum 
with n ~ 1 has been modified in-between, during the 
radiation-dominated era. Indeed, during this stage, den- 
sity fluctuations on scales larger than the Horizon keep 
growing whereas they oscillate on small scales, due to 
the pressure associated with the coupling to the radiation 
component of the Universe (photons). This implies that 
fluctuations SlO^) are multiplied by a transfer function 
T{k) that decays as at high wavenumbers. Then, the 
"initial" density power spectrum Psg (fc) used to study the 
formation of large-scale structures in the late Universe is 
the primordial one, with n ~ 1, multiplied by T(fc)^. This 
yields a curved cold dark matter power spectrum, with a 
local slope n that runs from 1 at low k to —3 at high k. 
Thus, today at z = we have n ~ —2 on galactic scales 
and n ~ — 1 slightly above cluster scales ^24|, [2^. This 
corresponds to the range studied in this article. Then, 
power-law power spectra with —3 < n < 1 can model the 
dynamics on the ranges of interest for specific purposes. 
On the other hand, the simplifications associated with 
power-law power spectra, such as the scaling laws p9p - 
(HT|) and the self-similar evolution (^5)) seen below, can 
be used to check the accuracy of numerical algorithms 
and to shed light on the dynamics |27i] . 



C. Spherically symmetric statistics 

In this article, in order to take advantage of the statis- 
tical isotropy of the system, we focus on two spherically 
symmetric quantities, the overdensity, rjr, and the mean 
divergence, Or, within spherical cells of radius r, which 
we define as 



Vr 



and 



m(< r) pr 
PoV po 



1 + 5r with Sr 



dx 

v 



e.^^/^^^(x) = -l/^dxu(x).x. 



5(x), 
(20) 



(21) 



where we used Eq.Q. Here V and S are the volume 
and the surface of the (d— 1)— sphere of radius r, x the 
unit radial vector, and we multiplied the divergence 9 
by time t in the definition (j2ip to have a dimensionless 
quantity O^. The moments (O^) can be understood as 
dimensionless spherical velocity structure functions, the 
usual longitudinal velocity structure functions being de- 
fined as ([(u(x) — u(0)).x]P) for a given direction x and 
length |x|, while in (pij) we integrate over all directions. 
In one dimension, d = 1, up to a sign, O^ is simply the 



dimensionless velocity increment over the distance 2r, 

d=l: e., = -L[u{r)-u{~r)]. (22) 
2r 

In arbitrary dimension, —Or is the dimensionless velocity 
increment over distance 2r averaged over all directions 
about a given point. We investigate in this article the 
probability distributions V{r]r) and 7^ (Or) in the quasi- 
linear regime (i.e. at large scales or early times), and 
their tails in any regime. The system being homogeneous 
we can focus on the cell that is centered on the origin, 
and this gives in Fourier space 



(23) 



Sr / dk(5(k)W^(fcr), with 



W{kr) 



dx 

v 



= 22 r 



Jd (kr) 

2 

(fcr)5 



In the linear regime we obtain from Eq.fT^ 



Oir 



t 

v 



dx Uo(x).x. 



(24) 



(25) 



For the initial conditions (fT5|) the linear density contrast 
(5ir is Gaussian, of mean zero, ((5i,r) = 0; and covariance 
C5i,,(»'i,f2) with 



C5i,,(ri,r2) = ((Jirif^Lrs) 



dk k'^-^Ps^ {k)W{kri)W{kr2) 



oc 



(ri-Hr2)"+3 



n + 3 d+1 



d+l: 



4rir2 



(26) 
(27) 

(28) 



where the last relation ([28)1 only holds for —3<n<d—l 
(if n > d — 1 the integral (177)) diverges at high k and the 
correlation Csj^^{ri,r2) is a distribution, such as a Dirac 
distribution for {n = 0, d = 1}). 
Then, the linear variance reads as 



n < d-2 



(29) 



Indeed, we can note that the integral (|27p converges at 
/c — > for any n > — 3 but only converges at fc — > cx3 for 
n < d — 2, when ri — r2. Therefore, in dimensions d < 3, 
the variance cr|^ shows an ultraviolet divergence for a 
power index in the range d — 2 < n < 1 (we only con- 
sider the range — 3 < n < 1 in this article). Of course, as 
soon as t > 0, the nonlinear evolution associated with the 
Burgers dynamics (shocks) makes the nonlinear variance 
(i5r) finite. In such a case one could also study the den- 
sity field smoothed by a Gaussian window, c>c e"'''' ^^^^ \ 
instead of the spherical top-hat ([^0]) . to obtain a finite 
linear variance crg^^, but we shall not investigate this al- 
ternative in this paper. 

We also introduce the spherical component of the ini- 
tial radial velocity, Mor, which reads from Eq. (|^5|) as 



Uor 



td 



(30) 



5 



n 


d 


C„o,,(ri,r2) 


C^0,.(»"l.»"2) 







1 


- r2) 


ri 


t'Soiri - r2)/(rir2) 





3 




ri(3r2 -ri)/(2r2) 




-1 


2 


ri/r-z 


rf[l +ln(r2/ri)]/2 




-2 


1 




r?('3r9 — ril/6 

' 1 V J- / / 




-2 


3 


ri-r?/(5ri) 


rf(r-?-5rir2 + 10r-i)/(20r2) 


9f2[l/r2-r?/(5ri)] 


n 


d ^ oo 


rir2(r? + ri)-("+^)/Vrf' 


[(r? + ri)(i-")/2 - r\- " r^"]/(d^(n2 - 1)) 





TABLE I: The initial velocity and potential correlations Cua^(ri,r2) and C^o^(ri,r2) of the spherical component of the initial 
conditions, for some values of n and d where Eas. (l27|l and (I33p simplify (with a dimensional normalization factor set to unity). 
The last column shows the covariance of the linear density contrast Slt within radius r at time t. Here we assumed ri < r2, 
except for the velocity and density correlations in the case {n = 0, d = 1}, and formulae for ri > r2 are obtained by exchanging 
ri and r2. The correlations are singular at ri = r2, except in the limit d ^ oo shown in the last line (where the velocity and 
potential correlations also become vanishingly small as compared with the density correlation). 



since V — {r/d)S. This is the mean initial radial velocity 
at radius r. As with rjr and O^, for spherical components 
we note the dependent coordinate r as an index, to dis- 
tinguish from the d-dimensional field Uq (x) (but contrary 
to rir and 0^, Ur is the mean radial velocity at radius r, 
rather than within the volume V). It will also be useful 
to consider the spherical component of the initial velocity 
potential, which we define from Eq.([3]) as 



dr' 



1 



Ar'r'SLr', (31) 



that is, we choose to normalize the initial potential by 
■00(0) — 0. Then, the initial radial velocity and potential 
two-point correlations are 



92 



rir2 



which can be obtained from Eqs.([26 |l -([28 |) . and 

2^+Vir(f + 1)^ 

/ J|_i(/cri) 
V (fcri)*-i " 

with a variance 

- 3 < n < 1 



(32) 



r(rf/2) cP J, 
2i-f \ / J±_i{kr2) 





'■4'Or 



, (33) 



(34) 



Note that ti^^^ is finite and well-defined over the whole 
range — 3 < n < 1. We give in Table |T] the initial ra- 
dial velocity and potential correlations Cu„,, (ri,r2) and 
C^Q^(r'i, for a few low integer values of n and d where 
they take a simple form, as well as the limit d ^ 00. We 
also show the covariance Cs^^ (ri , r2) of the linear density 
contrast at time t within radius r, from Eg. ([77)1 . We can 
check that they satisfy Eq. . The formulae are written 



for Ti < r2 (except for the velocity and density correla- 
tions in the case {n — 0,d — 1}) with a dimensional 
normalization factor set to unity. For these power-law 
initial power spectra, the normalizations used in Table U 
can always be achieved by a rescaling of spatial coordi- 
nates. 



D. Hopf-Cole solution and self-similarity 

As is well known, the nonlinear Burgers equation 
((D) can be solved through the Hopf-Cole transforma- 
tion ^, 'b'], by making the change of variable ip{x,t) = 
2z^lnS(x, t). This yields the linear heat equation for 
S(x, t), which leads to the solution 



7/'(x,i) = 2i/ln 



dq 



exp 



Then, in the inviscid limit v 
method gives P, Q 



V'(x,i) 



max 
q 



■00 (q) 



V'o(q) _ |x-qP 

(35) 

0+, a steepest-descent 



(36) 



2t 



If there is no shock, the maximum in p6p is reached at a 
unique point q(x, t), which is the Lagrangian coordinate 
of the particle that is located at the Eulerian position x 
at time t Moreover, this particle has kept its initial 

velocity and we have 



u(x,i) = Uo(q) 



q(x,i) 



t 



(37) 



If we have a shock at position x there are several degen- 
erate solutions to p6p and the velocity is discontinuous 
(as seen from Eq. (j37p . as we move from one solution q_ 
to another one q+ when we go through x from one side 
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of the shock surface to the other side) while the density 
is infinite. The solution (j36p has a nice geometrical in- 
terpretation in terms of paraboloids [J, [jj. Thus, let us 
consider the family of upward paraboloids 7^x,c(q) cen- 
tered at X and of height c, with a curvature radius <, 

^x,c(q) = ^^+c. (38) 

Then, moving down Px,c(q) from c = -|-oo, where the 
paraboloid is everywhere well above the initial potential 
-00 (q) (this is possible for the initial conditions ^ since 
we have |'!/'o(q)| ~ which grows more slowly 

than at large distances), until it touches the surface 
defined by ipQ (q) , the abscissa q of this first-contact point 
is the Lagrangian coordinate q(x, t) . If first-contact oc- 
curs simultaneously at several points there is a shock at 
the Eulerian location x. One can build in this manner 
the inverse Lagrangian map x i— > q(x, t) . 

For the initial conditions ^ that we consider in this 
paper, the rescaled initial velocity potential tpo{Xq) has 
the same probability distribution as X'^^'"'^^^^o{q) for 
any A > 0, when we normaHze by Uo(0) = and 00(0) = 
0, as seen in Eq.©. Then, the explicit solution gives 
the scaling laws 

0(x,i) i^^(<^x,l), (39) 
u(x,i) 'i?' t^u(t^x,l), (40) 
q(x,t) q(i^x,l) . (41) 

For the spherical overdensity r/r and the spherical velocity 
increment Qr this yields 

(42) 

that is, the distributions V{rjr\ t) and 7^(0^; t) of the over- 
density and velocity divergence at scale r and time t only 
depend on the ratio r/i^/("+'^^. 

These scalings mean that the dynamics is self-similar: 
a rescaling of time is statistically equivalent to a rescaling 
of distances, as 

A>0: t^Xt, x^A^x. (43) 

Thus, the system displays a hierarchical evolution as in- 
creasingly larger scales turn nonlinear. More precisely, 
since in the inviscid limit there is no preferred scale for 
the power-law initial conditions the only character- 
istic scale at a given time t is the so-called integral scale 
of turbulence, L{t), which is generated by the Burgers 
dynamics and grows with time as in (I43p . 

L(i) cxt2/(»+3). (44) 

It measures the typical distance between shocks, and it 
separates the large-scale quasi-linear regime, where the 
energy spectrum and the density power spectrum keep 



their initial power-law forms (fTTj) and (fT8|) . Ps{k,t) cx 
^2f^n+3-d^ from the small-scale nonlinear regime, which 
is governed by shocks and point-like masses, where the 
density power spectrum reaches the universal white-noise 
behavior (i.e. Ps{k,t) has a finite limit for k ^ 1/L{t)). 

This self-similar evolution only holds for n < 1 , so that 
|0o(q)| grows at larger scales, see for instance Eq.lfM]). 
and n > —3, so that |0o(q)| grows more slowly than 
and the solution ((36|) is well-defined [l3|. This is the 
range that we consider in this paper. The persistence of 
the initial power law at low k for the energy spectrum, 
E{k, t) cx fc"+^~'', that holds in such cases, is also called 
the principle of permanence of large eddies [l3l | . 



III. PERTURBATIVE EXPANSION AND 
ZELDOVICH DYNAMICS 

Although the Burgers dynamics can be integrated 
through the Hopf-Cole solution ([36]) , the computation of 
its statistical properties for random initial conditions re- 
mains a difficult problem for general n and d. Onl y in the 
peculiar one-dimensional cases n = [ij, [l^, [28l. l29l [30j 
and n = -2 ^14, 1^ [13, d, with d = 1, where the ini- 
tial velocity is a white noise or a Brownian motion, one 
can derive explicit analytical results, taking advantage of 
the Markovian character of these two specific stochastic 
processes. For general n and d one must resort to ap- 
proximation methods, such as perturbative expansions, 
as with most nonlinear dynamics. In particular, at early 
times we may look for the solution of the equations of 
motion ([l])-(l2]) as an expansion over powers of time. 




xGp(ki,..,kp)^o(ki)...^o(kp). (46) 

The Dirac factors express the invariance through trans- 
lations of the equations of motion, Fi = Gi = 1 from 
Eqs. pT)) . and the higher-order kernels Fp and Gp obey 
a recursion relation that is obtained by substituting the 
expansion ([45l) - (|46l) into the equations of motion ((T])-([2]). 
This yields in Fourier space 

pFp(ki, ..,kp) - Gp(ki, ..,kp) = 

E -M) Fp-^i^i+i, ..,k,), (47) 
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and 



(p~l)Gp(ki,..,k 



p-i 

E 



^ 2|ki,,|2|k,+i,p|2 

X (ki , . . , kf ) Gp_£ 1 , . . , kp ) , (48) 



where we note = ki + 1 
gives for p = 2 the kernels 



kj with j > i. 



i^|(ki,k2 



and 



G^(ki,k2 



(ki2.ki)(ki2.k2 



fc;2(ki.k2) 



This 



(49) 



(50) 



where we defined F|(ki,k2) = [i^2(ki, k2)+F2(k2, ki)]/2 
and G| the symmetrized kernels. In Eqs. (|45|) -(|50 |) we 
took the inviscid limit = 0+. Then, the effects of the 
infinitesimal viscosity (i.e. the formation of shocks) com- 
pletely disappear in these perturbative expansions. This 
implies that taking shocks into account requires non- 
perturbative methods. 

Note that the expansions (|45)) - ([50l) over powers of time 
are also expansions over powers of the initial velocity fluc- 
tuations 6*0 (x), or equivalently over powers of the linear 
density mode (5i(x, t) given in Ea. (|17p . Since the ampli- 
tude of the linear density fluctuations decreases at large 
scales, as seen in Eg.ip^. the perturbative expansions 
apply to both limits of early time or large scale. In par- 
ticular, in these limits the distributions V{rjr) and V{Qr) 
converge to the Gaussian of variance cr^^^ . Then, at early 
times, or when the system is smoothed over large scales, 
the displacements of particles are small and one recovers 
at leading order the linear theory of section Hi Bl that is 
set by the initial conditions. A simple example is pro- 
vided by the case {n = —2, d = 1} of one-dimensional 
initial Brownian velocity [T3, [l8'|. However, this only 
holds for n < d — 2, where the linear density variance 
as^^ is well deflned. For n > c? — 2 it is not possible to 
neglect shocks as soon as i 7^ 0, and the distributions 
V{rir) and V{Qr) remain far from Gaussian at any time 
and scale. This is for instance the behavior obtained in 
the case {n = 0, d = 1| of one-dimensional initial white- 
noise velocity [H, [H, [2^, [s^l • We shall recover these two 
different behaviors in the following sections. 

The perturbative approach (|45)) - ((46|) is the standard 
method used in cosmology to study the gravitational dy- 
namics at large scales and early times [U, [2^ |3l| (in 
this case the equation of motion ([T]) gets two new linear 
terms, associated with the gravitational force and a fric- 
tion term that comes from the expansion of the Universe 
and the change to comoving coordinates, but the non- 
linearity is the same and the perturbative expansion is 
similar). Indeed, in the standard cold dark matter model 
[2^ . the amplitude of the linear density fluctuations de- 
creases at larger scales (i.e —3 < n < 1 as in the present 
paper, with the same deflnition of the power-spectrum 



index n for d = 3), and the perturbative approach allows 
to describe the large scale structure of the Universe (e.g. 
beyond the scale associated with clusters of galaxies to- 
day), that is, the cosmic web formed by voids, filaments 
and walls that join the nonlinear high-density objects 
such as galaxies or clusters of galaxies. In the hydro- 
dynamical context, perturbative expansions over powers 
of time, such as (H5|) - (H5|) . have been used for instance 
in [H, [111 to study Eulerian and Lagrangian two-point 
correlations. They can also serve as a basis for Pade ap- 
proximants that attempt to improve the convergence of 
the series [H, H^l • 

From the point of view of the perturbative expansions 
(|15)) - ([5D)) . the Burgers dynamics ([T]) becomes equivalent 
in the inviscid limit to the Zeldovich dynamics [25|, ob- 
tained by setting the right hand side in Eq.([T]) to zero. 
This describes the free motion of collisionless particles, 
that always keep their initial velocity Uq and can cross 
each other. In a Lagrangian framework, the trajectory of 
the particle of initial Lagrangian coordinate q = always 
reads as 



x(q,t) 



tuo(q). 



(51) 



as in Eg. ([57)1 that only held before shocks. Before orbit- 
crossings the conservation of matter gives for the density 
field 



p(x)dx = podq, hence l-|-5(x) 
This gives 



c^x 

dot ( — 
oq 



l + (5(x) = y"dq (5n[x-q-fuo(q)], 



(52) 



(53) 



which remains valid after orbit crossing as we integrate 
over all streams that pass through position x at time t. 
In Fourier space this yields 



dq 

{2ttY 



g ik.q ^g-ik.tuo(q) 



1 



(54) 



Then, expanding the exponential over Uq directly gives 
the s ym metric kernels Fp associated with the expansion 

65]) Hi, 



F;(ki,..,kp) = - 



1 ki^p.ki ki^p.kp 



kl 



(55) 



which agrees with (|49p for p — 2. From the perturba- 
tive expansions (|^ - ([5n)l we can compute the cumulants 
of the smoothed density contrast rjr and velocity diver- 
gence Qr (in the quasi-linear regime where shocks do not 
contribute, that is, leading-order terms at early times and 
large scales for n < d — 3, as discussed in section |IV] be- 
low). For instance, substituting the expression ([IS]) . the 
density three-point correlation reads in Fourier space as 



((S(ki),5(k2)5(k3))c = (5D(ki-Hk2+k3) 

X [2P5j/ci,t)P5j/c2,t)F|(ki,k2) +cyc. 



(56) 
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where "eye." stands for two terms assoeiated with eyehe 
permutations over {ki,k2,k3} of the previous term, 
while the dots stand for higher-order terms. Then, from 
Ea. ip^l the cumulant of order three of the overdensity 
within radius r writes 

{Vr)c - 6 J dkidk2P5jfcl,i)F5jfc2,i)^^|(kl,k2) 

xW{kir)W{k2r)W{\ki + kajr) + ... (57) 

Using the properties of Bessel functions, such as their 
addition theorem, one obtains for instance in dimension 
rf = 3 SHI, 

{rjf.), = {l-n)<jl^ + ... (58) 

One can use this method to derive the leading-order term 
of all cumulants {'rj^)c and {Qr)c- Then, from the charac- 
teristic function if{y), defined from the Taylor series (jM]) 
below, one can reconstruct the distributions V (vr) and 
■p(9r) in the quasi-linear regime, cts^^ <C 1 [1^, [s^. We 
shall describe in section IIVI below another method that 
directly gives the generating function (p{y) without using 
the expansions ((45l) - ((50)) and that allows to go beyond 
the singularities associated with the Taylor series ([M)) . 
Let us recall here that the previous results only hold for 
the case n < d — 2, where the linear theory is meaningful 
(i.e. (TSj^^ is well-defined). 

IV. QUASI-LINEAR LIMIT 

A. Distribution of the density within spherical cells 

We consider here the probability distribution, V{r]r), 
of the overdensity rjr within spherical cells introduced in 
(PO)) . More precisely, we investigate its asymptotic form 
in the quasi-linear limit, defined as cfsl,- ^ 0- Therefore, 
we restrict ourselves to the range n < d — 2 (in addition 
to — 3 < n < 1) so that the linear variance cr|^ is well 
defined, see Eg. ipS]) . Taking advantage of the statistical 
isotropy of the system, we apply to the Burgers dynamics 
the steepest-descent method (instanton technique) that 
was devised in [39j] for the coUisionless gravitational dy- 
namics. 

1. Action S[Sl] 

To obtain the quasi-linear limit of the probability dis- 
tribution V{rjr) it is convenient to first introduce the mo- 
ment generating function '^{y), 

*(y) - (e-^"") = / Ai^re-y^^Vivr), (59) 
Jo 

from which we can recover V{r]r) through the inverse 
Laplace transform 

'PM= / T^ey'^'-fiy). (60) 



Since the system is fully defined by the Gaussian linear 
density field at the time of interest, 5l(x) (we usually 
omit the time dependence as t can be seen as a mere 
parameter, since we only consider equal-time statistics), 
the average ([59)) can be written as the path-integral 

*(y) = (detC7/)i/2 I v5l e-^'^I^-l-^^-c-.-^.s,^ ^^^^ 

where 77r[^L] is the functional that affects to the initial 
condition defined by the linear density field <5l(x) the 
nonlinear overdensity 77^, built by the Burgers dynamics 
([I])-® at time t, within the spherical cell of radius r 
centered (for instance) on the origin x = 0. Here and in 
the following we use the short-hand notation for scalar 
products 

5L-Ci^.5L= y"dxidx2(5L(xi)Q-i(xi,X2)5L(x2), (62) 

where C^^ is the inverse of the two-point correlation (|19p . 
Equation ([61) is exact but the difficulty of the problem 
is hidden in the nonlinear functional r^rl^L]- In order to 
make some progress, we consider the quasi-linear limit, 
'^Slt ~^ 0: associated with large scales or early times. 
Then, it is convenient to rescale the moment generating 
function as [soj 

^-(y) = e-'^^^'^'^^^/'^'i^, (63) 

where (p{y) is the cumulant generating function, which 
has the Taylor expansion 

Substituting Eq. ([^ into Ea. ([CT|) gives 

e-v(?/)/<. = (det C^^f'^ j V5l g-^I^^'/^. , (65) 
with the action S[6l\ given by 

S[5l] = VvASl] + ^ 5l.CJ^.5l (66) 

The rescaling ([55]) allows us to derive the quasi-linear 
limit through a steepest-descent method, since the ac- 
tion S[5l\ no longer depends on the amplitude of the 
two-point correlation Cg^^ (as cx C^^ ) and the path- 
integral ()65|) is clearly dominated by the minimum of the 
action S in the limit ag^^^ 0. 

Here we may note that the use of a path-integral for- 
malism to analyze dynamical systems such as Eqs.(|T|)-([2]) 
is a standard approach, following the operator formalism 
of Martin-Sig gia- Rose px| or the functional method of 
Phythian j4l[. |42|. l43j . In such a framework, the path- 
integral (|6ip is rewritten in terms of the nonlinear fields 
(5(x, t) and 0(x, t), so that rjr is a simple linear functional 
of S as in ([^D|) , by introducing a Dirac functional such as 
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SoidtU + (u.V)u — i/Au — (and similarly for 6) to en- 
force the equation of motion ([1]), where, depending on the 
system, ^(x, can represent both a stochastic external 
forcing and the random initial conditions. Taking care of 
the Jacobian, which is usually equal to unity thanks to 
causality 44], one obtains a path integral such as ([65|) . 
but over the nonlinear density field (5(x, t) and its con- 
jugate A(x, i), rather than over 6l{x.), and over the ve- 
locity pair |0(x t), jn(x, t)}. This procedure is described 
in details in f3J| for the Zeldovich dynamics recalled in 
section Hn] above, that amounts to set the right hand side 
in the Burgers equation ([1]) to zero, see also jl^ for the 
collisionless gravitational dynamics. For noisy dynamics, 
where one adds a stochastic external forcing, this method 
is presented for instance in [i^ |43| for the forced Burgers 
dynamics and in ^] for the forced Navier-Stokes dynam- 
ics. 

Then, one obtains a cubic action S[5, 6; A, /i]. Expand- 
ing over the cubic term gives back the perturbative re- 
sults discussed in section lllll as one recovers an expansion 
over powers of the initial power spectrum Pgj^ . On the 
other hand, this path integral can serve as a basis for 
other expansion schemes, such as large- methods [13], 
that recover at leading order Kraichnan's direct interac- 
tion approximation when applied to the Navier-Stokes 
equations [igllsol]. 

Alternative expansion schemes, where one does not ex- 
pand over powers of some coupling constant or parame- 
ter, are provided by steepest-descent methods (instanton 
techniques [3]) where one expands around a saddle-point 
of the action S. If this saddle-point is non-perturbative 
this approach can go beyond perturbative expansions 
such as those described in section IIIH as we shall see 
more clearly in section |V] below. This approach has been 
applied to the forced Burgers dynamics in [i^, \^ and 
to the forced Navier-Stokes dynamics in [isj . In particu- 
lar, this allows to obtain the right exponential tail of the 
probability distribution of the velocity increment [46l . Isil ] 
and its left power-law tail [s^ . 

In the noiseless case, a problem that arises when 
one tries to apply this method to the standard action 
S[S, 9; A, fj], directly obtained from the equations of mo- 
tion as described above, is that this action is highly sin- 
gular when there is no external forcing. Indeed, in such 
cases the dynamics is fully deterministic, so that the sys- 
tem is fully defined by the initial condition 0o(x) (or by 
(5l(x) at a given reference time). Then, one can check 
that this action S is only finite for fields that obey the 
equations of motion (Il])-((2]) and infinite elsewhere, which 
simply means that the path integral only counts fields 
that are solutions of the dynamics, as it should. There- 
fore, the action is only finite over a lower dimensional 
subspace parameterized by ^o(x), that is, the time de- 
gree of freedom of the fields (5(x, t) and 0(x, t) is not real. 
Then, the action has no finite second-derivative and the 
steepest-descent approach is not very well defined. More- 
over, any expansion point must be an exact solution of 
the dynamics so that this approach does not bring much 



progress. 

By contrast, the path integral (|65p only involves the 
true degrees of freedom of the system, parameterized by 
the linear density field <5l(x, t) at the given time of inter- 
est (i.e. we do not integrate over non-existent time de- 
grees of freedom). Then, the action S[6l] is finite and has 
a well-defined second derivative, at least close to y = 
and 6l = 0, so that the steepest-descent approach rests 
on firm grounds. Moreover, the difficulty associated with 
the nonlinear functional i]r[SL] would not be overcome by 
using the standard action S[S,9; X, fi], since in this case 
too we would need to study exact solutions of the dy- 
namics. Thus, the action ([66|) is well-suited to the appli- 
cation of the steepest-descent approach to deterministic 
dynamics, as we shall see in the following. In practice, 
in order to handle the term r?r[^L], one must be able to 
obtain saddle-points where the dynamics can be explic- 
itly solved in simple terms. In our case, it is natural to 
take advantage of the statistical isotropy of the system to 
look for spherically symmetric solutions of the dynamics. 
Then, this requires to focus on spherically symmetric ob- 
servables, such as rjr and 9^ defined in Eas. (|20l) - ((2T|l . so 
that spherical initial conditions can also be saddle-points 
of the action 5J(5l]. In fact, in such a case, a minimum of 
the action with respect to spherically symmetric initial 
conditions is automatically a saddle-point with respect 
to non-spherically symmetric initial conditions (but not 
necessarily a minimum). In our case, we shall see below 
in section llV A 2l that we really obtain a local minimum in 
the quasi-linear regime (i.e. for small y and Sl). Then, 
even if there exists another local minimum reached for 
some non-spherical initial conditions, which requires a fi- 
nite fluctuation Sl, such a contribution is exponentially 
subdominant in the quasi-linear limit, ct^j^,, — > 0, so that 
we obtain exact results in this limit, without the need 
to explicitly study the functional r^rif^L] over all possible 
non-spherical states. 

Note that for non-spherically symmetric observables 
(for instance we could choose cubic cells to define the 
mean density rjr and velocity divergence 0^) we could 
apply the same approach and look for local minima with 
respect to spherically symmetric initial conditions. How- 
ever, these would no longer be saddle-points with respect 
to non-spherical fluctuations, and would not be the true 
minima of the action. Then, one would only obtain lower 
bounds for the asymptotic behavior of the distribution 
ViX), where X is any observable and is not necessarily 
spherically symmetric, in the quasi-linear or rare-event 
limits. One can expect that the exponents obtained in 
this fashion would remain correct, but the numerical fac- 
tors within exponential tails would only be approximate. 
The asymptotic behaviors obtained within this approxi- 
mation would clearly show the same qualitative proper- 
ties as those obtained for the spherical observables stud- 
ied here (since one uses the same initial states). Thus, 
it is straightforward to obtain lower bounds for any dis- 
tribution V{X), in the quasi-linear or rare-event limits, 
from the method described in this article and we restrict 
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ourselves to the spherically symmetric observables ((20 
(PT|) in the following. 



2. Spherical saddle-point 

As explained above, and as for the gravitational dy- 
namics [39], taking advantage of the spherical symme- 
try of the action we can look for a spherical 
saddle-point. Indeed, since the first functional deriva- 
tive, VS /VSLiti), taken at a spherical linear density field 
(5L(q), is spherically symmetric, it only depends on |q|. 
Then, the first variation AS due to a non-radial pertur- 
bation A(5L(q) vanishes. 



AS= / dq 



ASUci) = 0, 



when 



SUci) ^ SlIM) and / dqA,SL(q)=0, 

J\q\=q 



(67) 



(68) 



where the second equality is the integral over angular 
variables at any radius q. Therefore, a saddle-point with 
respect to spherically symmetric states (i.e. radial de- 
grees of freedom) is automatically a saddle-point with 
respect to angular degrees of freedom, whence a true 
saddle-point with respect to any infinitesimal perturba- 
tion A6l{<i)- Then, we can restrict the action S[Sl] to 
spherically symmetric initial conditions and look for its 
minimum within this subspace. For such initial condi- 
tions, the action can be expressed in terms of the one- 
dimensional field SLq', defined as in Eq. ipD]) over < 
q' < cxD, (we note by the letter q initial Lagrangian radii, 
to distinguish them from the Eulerian radii r reached at 
time t). This reads as 



S[SLq'] ^ yVr[SLq'] + 



■^S,q'^.Ci^],.S,,^, (69) 



where q' is a dummy variable and Csj^^{q'i, q2) is the co- 
variance introduced in Eq. (j26p . Then, saddle-points of 
the action are given by the condition VS/VSLqi = 
over < g' < cxd, that is. 



VS 



Lq' 



dq"CiMq',q")dLq" =0. (70) 



Multiplying by the operator Ca^^^ this reads as 

SLq' = 



dq Cs^M,1 • (71) 



'VSLq" 



Next, we note that if there have been no collisions (i.e. no 
shocks) until time t, the spherical collapse or expansion 
has remained well ordered, and the mass m within the 
radius r comes from the matter that was initially located 
within a Lagrangian radius q at time t = 0. Then, the 
overdensity rjr = m/{pQV) is also given by r/r — (q/r)'^ 



and it only depends on the initial Lagrangian coordinate 
q of the shell that is located at radius r at time t. On the 
other hand, in the inviscid limit the Burgers dynamics ([1]) 
implies that particles that have not collided yet have kept 
their initial velocity Uq. Therefore, for a spherical state 
the initial Lagrangian radius q is related to the Eulerian 
radius r hy r = q + t uoq , whence 



Vr = [q/rf = 1 + 



with 



nsLq) = ( 1 - % 



= T{5Lq). (72) 



(73) 



where uoq is the initial radial velocity at radius q and 
we used Eg. ([50)1 (for spherical initial conditions we have 
Uo(x) — upajx). Thus, the overdensity rjr only depends 
on the initial velocity at the Lagrangian coordinate q, 
whence on the linear density contrast SLq within the La- 
grangian radius q. As a consequence, it is independent of 
infinitesimal perturbations to the initial profile 6Lq' over 
inner or outer shells {q' < q or q' > q), that only redis- 
tribute matter at smaller or larger radii. On the other 
hand, under an infinitesimal perturbation AS^q' the La- 
grangian radius q and the overdensity rjr are modified as 
q ^ q + Aq and rjr ~^ rjr + A7]r . From Eq. ([7^ we obtain 
at first order, 



Anr = J"{5Lq) 



d5 



Lq' 



dq' 



Aq + ASLq 



Ar]r 



Aq 



(74) 
(75) 



This leads to Arjr oc AS^q, which means that the func- 
tional differential 'Drjr/'DSLqii in Ea. (|7T|) is a Dirac distri- 
bution centered on q" = q, in agreement with the previ- 
ous discussion, and we directly obtain the initial profile 
of the saddle-point as 



SLq' oc CsL^iq',q) whence SLq' = SLq 



^1 

OLq 



(76) 

Using Eq. ([50)1 , this also gives for the initial velocity pro- 
file 



UQq' 



q SLq' 
UOq - -J— 
q OLq 



UOq ■ 



CuoM'^q) 



(77) 



Next, the amplitude SLq, or the Lagrangian coordinate 
q, can be determined by substituting the profile ([76)) into 
the action ()69p and looking for its minimum with respect 
to SLq- This reads as 



S = vT{SLq) 



2a2 



(78) 
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Then, defining the variable t and the function Q (r) by 



= -<5 



Lq- 



G{t) = T{5Lq) = 77r 



the action ([78|) and its derivative read as 
-2 dS 



S = yg(r) 



dr 



vQ' 



(79) 



(80) 



Therefore, since at leading order in the quasi-linear limit, 
the cumulant generating function ^p{y) is given by the 
minimum of the action S[5l\ from Ea. (|65p . it is given by 
the implicit system 



v{y) = yGir) + — with T 



(81) 



Thus, the generating function (/^(y) is also the Legendre 
transform of the function — r(t/)2/2, as defined by 



(^(y) = min 



mm 
Q 



vQ 



rig)' 



(82) 



To make sure that the solution ((8T|) is indeed relevant, 
we must check that it is indeed a local minimum of the 
action (and not a maximum), in agreement with (j82p and 
the original path integral (|65p . This directly follows from 
the expression (|66|) . Indeed, for y = the saddle-point 
obtained above is simply Sl = 0, i.e. r = 0, and the 
Hessian of the action at this point is a^^ C^^ which is 
strictly positive. Then, by continuity, for small y the Hes- 
sian around the saddle-point given by Eq. (|ST|) is positive 
which ensures that it is a local minimum. As we shall see 
below in section HV A 31 for some cases it may only be a 
local minimum, but the global minimum associated with 
finite density contrasts is irrelevant in the quasi-linear 
limit: it corresponds to the tail of the distribution V{rjr) 
and it is exponentially suppressed in the limit cr^^^, 0. 

It is clear that the procedure described above must re- 
cover the results that would be obtained for the leading- 
order term of the cumulants (??J?)c (which is of order 
a^g^ that ip{y) has indeed a finite quasi-linear limit 

in Ea. (j64p ) from the perturbative expansion presented 
in section Hill Indeed, in both cases we obtain an expan- 
sion over powers of tr^^^ (in the steepest-descent approach 
subleading terms would be obtained from Eq. (|65p by ex- 
panding the action around its saddle-point and perform- 
ing the Gaussian integrations), as we actually start from 
the unperturbed solution 5l — In fact, as shown for 
the case of the three-dimensional gravitational dynamics 
[53| . it is possible to derive the quasi-linear generating 
function Lp{y) from the perturbative expansion (j45p - (|46p . 
using its Taylor expansion (j64p and the leading-order 
term of each cumulant (?7^)c- This gives back Lpiy) as the 
solution of the implicit system |(5T|) for the unsmoothed 
case, where one has Q{t) = P{t) (so that there is no de- 
pendence on the initial conditions). Then, one can show 
that the same result is obtained in Lagrangian space, 
where the perturbative expansions are built in terms of 




q'/q 



FIG. 1: (Color online) The linear density profile of the spher- 
ical saddle-point for the cases n = —2, and 0.5, in dimension 
d = 3. This shows the integrated linear density contrast 5^^/ 
within the sphere of radius q' , from Ea. (|83|l . and q is the ini- 
tial Lagrangian radius of the shell that is at the radius of 
interest r at time t. 



the Lagrangian coordinates, and the shift from J-{Sl) to 
^/(t) as in Eq. (|79p is obtained through a mapping from 
Lagrangian to Eulerian space [13] . 

As shown in [soj and described above, the steepest- 
descent approach provides the quasi-linear generating 
function in a more direct fashion. In particular, the inte- 
grations over angles, as in Eq. ([57)) . are automatically in- 
cluded in the spherical dynamics, and the dependence on 
n, associated with the mapping from Lagrangian to Eule- 
rian coordinates (through the ratio crs^^/ag^^ in Eq. lfTS]) 
that measures the ratio of initial power on scales r and q) 
is automatically provided by the form of the action 5[(5l]. 
In addition, the definition of f{y) as the Laplace trans- 
form ()65p allows to give a meaning to ip{y) beyond the 
radius of convergence of its Taylor series. Finally, as dis- 
cussed below, the profile ((76l) of the saddle-point allows 
us to check the range of overdensity r]r and Laplace con- 
jugate y to which these results apply. Indeed, they only 
hold as long as the saddle-point has not formed shocks 
yet, which can only be checked from the knowledge of 
Eq.®. 

For the power-law power spectra (|18p the radial linear 
profile ([7S)l of the saddle-point reads as 



1 



2Fl(2±2,d±l;d+l;l) 



, (83) 



with X — q'/q, see Eg. ((28)) . This also gives the linear 
velocity profile through Eq. ljSOl) . We show in Figs. (TJ 
[21 the density and velocity profiles obtained in the three- 
dimensional case for n = —2, and 0.5. For integer values 
of n and odd d+n the hypergeometric function in Eq. (|83p 
simplifies and we give in Table [TTl the explicit forms ob- 
tained for some low dimensional cases. Note that the 
profile is singular at the Lagrangian radius q, associated 
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qVq 



n 


r?r > 1, e,. > 


r]r <1, Or <0 


n > d — 3 


shock as soon as t ^ 


-2 < 71 < d - 3 


no shock 


shock below a threshold 


-3 < n < -2 


no shock 


no shock 



TABLE III: This Table shows whether the saddle-point jTe]) 
forms a shock after a finite time, which corresponds to a finite 
threshold for the density r/r or the velocity divergence Or- We 
only consider the range —3 < n < 1 and d>l. If n>ci — 3 
shocks form as soon as t 7^ so that the saddle-point (|76p is 
never valid (but it should give a reasonable approximation if 
n < d - 2). 



FIG. 2: (Color online) The linear velocity profile of the spher- 
ical saddle-point for the cases n — —2, and 0.5, in dimension 
d = 3, as in Fig.[l] This shows the initial radial velocity uo{q') 
at Lagrangian radius q' . 
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1 

X 
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1 
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-2 


3 




5x^-1 


x{5-x^) 


5x^-1 


4 


4x^ 


4 


Ax'-' 


n 


00 


(1^ 


^J-(n+3)/2 







TABLE II: The linear integrated-density and velocity profiles 
of the spherical saddle-point for some values of n and d, where 
Ea. (|83|l simplifies. Here x = q' /q, where q is the Lagrangian 
radius associated with the Eulerian radius of interest r. The 
profile is singular at 2: = 1, except in the limit of infinite 
dimension, d — » 00, shown in the last row. 



with the radius r of the Eulerian cell. We also give in the 
last row the simple profile obtained in the limit d 00, 
where the singularity disappears. Note that in this infi- 
nite dimensional limit the profile still depends on n, in 
agreement with the fact that the profile decays as 
^-(n+3) large distance, independently of d. 

As expected, for all values of n in the range — 3 < n < 1 
that we consider in this paper, the density contrast van- 
ishes at large distance. For n < d — 3 it is monotonically 
decreasing but for n > d — 3, which corresponds to sig- 
nificant initial power at high wavcnumbers, it shows a 
peak at radius q. On the other hand, the radial velocity 
vanishes at the center q' — 0, in agreement with spher- 
ical symmetry, but it only decays at large distance for 
n > — 2. For n < — 2 it keeps growing at large distance 
(note that the initial velocity field only shows homoge- 
neous increments for n < — 1, so that this growth is not 



surprising) . 

To make sure that the saddle-point obtained above 
is relevant we must check that no shocks have formed, 
so that Eg. ((72)) is valid. The naive Lagrangian map, 
X = q + tuo(q), shows that a shock occurs when 
det(9x/(9q) — 0, that is when 1 -I- tduoq'/dq' = for the 
spherical saddle-point. The profiles show a singularity at 
radius q of the form \q' — q\'^^'"-^'^ . Then, for n > d — 3 the 
velocity has a spike at radius q with infinite left and right 
derivatives, so that shocks appear as soon as f 7^ 0. For 
n < d— 3, \uoqi I shows a sublinear (or linear) growth with 
q' hence there will be no shock, except at the center, for 
overdense saddle-points (particles reach the center before 
duoq'/dq' reaches —l/t). For underdense saddle-points, 
a shock appears after a finite time for —2 < n < d — 3, 
while no shocks form for —3 < n < —2 since the radial 
velocity grows with radius. We summarize in Table IIIII 
these behaviors associated with different ranges of the 
index n of the initial energy spectrum. 

Since in the quasi-linear limit we only probe small den- 
sity fluctuations we can use the saddle-point obtained 
above for n < d — 3, as shocks only appear after some 
finite time (or never). For n > d — 3 we should modify 
the saddle-point to take into account shocks. However, 
since for moderate times and density fluctuations this 
should only change the profile close to the Lagrangian 
radius q and give small modifications to the quasi-linear 
generating function (p{y) we shall keep Eq. (|8T|) below for 
n < d—2 (for n > d—2 where the linear variance (75^^^, di- 
verges it is not possible to neglect shocks) . Note that for 
large d this problem disappears, see also the last row of 
Table HIl For all cases shown in Table [U the saddle-point 
obtained above is relevant in the quasi-linear regime, as 
can be checked from the explicit forms of the velocity 
profiles and in agreement with the previous discussion. 
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3. Cumulant generating function tp(y) 

Close to the origin y — the solution of the implicit 
system (|8ip always satisfies 



0, y^O: g ^l-T, y 



T and (/s-^y-y. 



(84) 

Keeping only these low-order terms corresponds to the 
linear regime and gives back the Gaussian of variance 
(t|^ for the probability distribution 7^(?yr) on very large 
scales and early times. This agrees with the fact that the 
initial conditions are Gaussian and it means that over 
large scales or at early times, that is when (t|^ <C 1 
and for \6r\ ^ 1, we recover linear theory and the prob- 
ability distributions are still governed by the Gaussian 
initial conditions. Of course, this breaks down for the 
cases d — 2 < n < 1 where the linear variance cr^^ itself 
is not well defined, so that even on large scales the prob- 
ability distributions are strongly non-Gaussian, as seen 
for instance in [H [H, [H, [p| for the case {n = 0, d = 1}. 

For the case n < d — 3, where the quasi-linear regime 
considered in this section applies, the deviations from 
the Gaussian, associated with higher-order cumulants, 
appear for finite density contrast 5r, that is for finite y 
and T. Note that this corresponds to rare events, as ex- 
pected for a steepest-descent approach to be valid. Thus, 
in the quasi-linear limit we are sensitive to small but fi- 
nite values of {5^, y, t} around zero. Moreover, at leading 
order the moments and cumulants of the overdensity are 
given by the expansion ([M)) of (p{y) around the origin. 
In particular, the solution ([5T|) directly gives the leading 
order value of {ri^)c, that can also be derived from the 
perturbative expansion of the equations of motion 
described in section [Till 

The previous results are valid for any initial energy 
spectrum, provided asLr ''^sll defined. For the power- 
law power spectra ([T5| . using Eq. p5)) . we obtain from 
the definition ([75)l . 



3L9 



whence 



1 + -^-^ 

d 



g{T)^T[-Tg-^ 

In terms of the inverse hmction t(CJ) this reads as 

T{g)^d(g^ -g^ 



(85) 



(86) 



(87) 



so that is given by the parametric representation 



^ = 2 



{d-n-\)g^ - 2(d-n-2)g- 
+ (d-n-3)g'^l , 



(n+\)g- 



-2)g 



+ (^+3)^"^ 



(89) 

Expanding around t — 0, y ~ and g — 1, we obtain 
the series expansion of >f{y). Comparing with Eq. (j64p 
we obtain for the third and fourth-order cumulants in 
the quasi- linear limit. 



= 3- 



d-n-2 



(90) 



which agrees with ((58| for d = 3, and 



54 = 



in 



2\3 

Tic 



83 -f led^ ^- 84n ^ 21n2 - 36d(n -I- 2) 
d2 ■ 



(91) 

Note that although ^3 and ^4 as defined above are called 
the skewness and the kurtosis in the cosmological litter- 
ature, they are not exactly the skewness and the kurtosis 
defined in standard probability theory, the latter being 

defined as (?7r)c/('7r)c^^ and (?7r)c/('?r)c- The reason for 
the use of ([101)- (EJ) is that these quantities have a finite 
value in the quasi-linear limit discussed above (as seen in 
,3ft] they also have a finite value in the small-scale limit, 
associated with the highly nonlinear regime). 

For small integer values of n and d we can obtain ex- 
plicit expressions for the solution of the implicit system 
([88ll -([89 l) . by solving for g(ij) and substituting into <p(5). 
We give in Table HVl our results for a few such cases. Note 
that from the meaning of the variable g as the overden- 
sity rir within the radius r for the spherical saddle-point, 
see Eq.([79l), '^{y) is a priori determined by Eqs. ([88l) - ([89l) 
by letting g vary over the range < < 00. For the 
cases {n = — l,d = 2} and {n = — 2,d = 1} we obtain 
a generating function (^(y) that shows a singularity ys 
on the negative real axis, with y^ ~ —2 or —1/2, and 
the range y > ys corresponds to the full range g > 0. 
In the complex plane, there is usually a branch cut for 
y < ys or a pole at y^. For the cases {n = 0, d = 3} 
and {n — — 2,d = 3} it happens that the function y{g) 
is no longer monotonic over < ^ < 00 so that the in- 
verse g{y) is bivaluate and f{y) shows two branches. We 
show in columns 3-5 of Table HV] the quasi- linear branch, 
that contains the point {y = 0,g — l,ip = 0} and cor- 
responds to moderate density fluctuations. Columns 6-8 
show the second branch that corresponds to large fluctua- 
tions (very low densities, < g < 1/8, for {n = 0, d = 3}; 
very high densities, 2^/2 < g < 00, for {n — —2, d = 3}). 

For the general case, the behavior of the cumulant gen- 
erating function (p{y) defined by the implicit system ([8T1) 
and the presence of singularities can be obtained from the 
asymptotic behaviors at large and small overdensities g, 
using Eqs. ([87l) - ([89l) . For large densities we obtain 



-00 : r ~ ~dg ^ y , 
d(d — n — 3) 

V 7, 



d{n + 3) ,i+3^d 
7, Q I 



'2y 



_d(n + 3) 



(92) 
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27-9V9-6!/ 


9 ^ ^3 


i < < OO 


27+9V9-6!/ 
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-1 
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1+I//2 


— 2 < y < oo 


< e < OO 








-2 
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Vl + 2y-l 


< y < oo 


< 5 < oo 








-2 


3 




-| < y < oo 


< 5 < 2\/2 




-| <y<o 


2\/2 < e < oo 






^3+^9+24:/ 


^3-^9+24;/ 


n 


oo 


(/3 = T+^, y = re^, C7 = e^; -■i<y<oo, 0<e<e; -\ < y < Q, e < G < oo 



TABLE IV: The cumulant generating function ip(y) of the overdensity 77^, in the quasi-linear Umit crs^^ — > 0, for a few values of 
n and d, where exphcit solutions of the system (|88p - (l89| l can be obtained. The columns 3-5 show the quasi- linear branch, and 
the associated range of {y,G} that contains the point {0, 1}. The columns 6-8 show the second branch, associated with very 
rare events, that appears in some cases. The last row shows the infinite-dimensional limit, d — > 00, which no longer depends 
on n but has no explicit form and shows two branches. 



Thus, we have two possible behaviors for Q +00, 
n > d — 3 : t —00, y —00, ip — )■ —00, (93) 



<d-3 



-00, y 



^ - 



+00. 



(94) 



As explained in section IIV A 21 and Table IIIIl the saddle- 
point approach studied here only exactly applies to n < 
d — 3 as shocks form for n > d — 3. However, we mention 
the case n > d— 3 in because this method should still 
provide a reasonable approximation for d—3 < n < d—2. 
The case n — d—S shows an intermediate behavior, as in 
the limit Q +00 we obtain y —d'^/2, and (p ^ —1 
ii d = 1 or (fi ^ —00 if d > 1. Thus it is closer to the 
case n > d — 3. Then, we can see that for n > d — 3 
larger densities are associated with more negative y and 
if and the function ipiy) is regular and monotonically 
increasing over ] — 00, 0] (or ] — d^/2,0]). This behavior 
is shown by the case {n = 0, d = 3} in Fig. [3] For 
n < d — 3, since from (|84p we have y — a.t Q ^ 1, the 
limit y ^ 0^ for large densities implies that the function 
y{G) is not monotonic over G € [l,+oo[ and shows a 
minimum ?/s < at some value t/s > 1. Around this 
point we have y — ys oc {G — Gs)^- This gives rise to 
a square-root singularity \/y ^ ys for the function (p{y), 
which shows two branches going from this point. A first 
branch goes through the point {y = 0, = 0}, it is the 
branch associated with moderate fluctuations, below Gs, 
that is most relevant in the quasi- linear limit. The second 
branch is associated with large overdensities above Gs- 
This behavior is shown by the case {n = —2, d = 3} in 
Fig.H 

For low densities we obtain 



: r ^ dt/ 2d ^ y ^ 
d{d - n - 1) 



d{n + 1) n+i-d 

G 

2 

-2?y 



d(n + 1)_ 



(95) 



Since we assumed n < d—2, so that the linear variance 
cr? is well defined, we have n-|-l — d< —1 and this 



gives rise to the two behaviors: 

n < —1 : T —>■ +00, y —>■ +00, ip — > -foo, (96) 



n > -1 



0+, y 



(97) 



Thus, for 71 < —1 the function (p{y) is regular and mono- 
tonically increasing over [0,-|-oo[ (case {n = — 2,d = 3} 
in Fig. [3]) while for n > — 1 it shows a square-root sin- 
gularity at some finite value ys > 0, associated with an 
undcrdcnsity Gs < ^, from which two branches leave (case 
{7i = 0,d = 3} in Fig. El . 

In the large-dimension limit, d — > 00, we obtain from 
Eqs.jTS]), (ESI), J^((5l) = e^^ and r -Sl- This gives 
the parametric representation of ip{y) shown in the last 
row of Table IIVI Since < G < 00 corresponds to 
—00 < T < 00 and ?/(r) has a minimum at Ts = — 1 
this generating function (p{y) shows two branches. The 
comparison with appendix A of [s^ shows that in the 
quasi-linear limit this leads to a log-normal distribution 
'P{rjr) (contrary to some statistical models, used to de- 
scribe fully developed turbulence, this is not related to 
some underlying multiplicative cascade process). Note 
that the dependence of </?(?/) and Vijir) on n disappears 
in this limit d — > 00, even though the density profile of 
the associated saddle-point keeps a dependence on n, see 
last row of Table [TTl 

As noticed in [18], it happens that in the case {n — 
— 2, d = 1} the quasi- linear result shown in the fourth row 
in Table |TV] actually gives the exact cumulant generating 
fimction defined as ^(y) — — ^Sp{—yY/p\, with Sp = 

{v?)c/{''lr)'c' ■ Note that for the quasi-linear limit we 
defined (p{y) in Eq. ((64| using (t|^_ instead of (?7r)c in 
the coefficients Sp, which is only equivalent at leading 
order. In this case {n — —2, d = 1} we actually have 
the exact equality {r]'^)c = '^SLr^ generic 
cases we expect the exact generating function ^(j/) and 
the variance (?7r)c to deviate at small scales or late times 
from their quasi-linear limits (p{y) and . 
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FIG. 3: (Color online) The density cumulant generating func- 
tion ip{y) in dimension d = 3 for the power indices n — 
(solid line) and n = —2 (dashed line), from rows 2 and 5 
of Table IIVI In both cases there is a singularity on the real 
axis, with i/s = 3/2 for n = and ys — —3/8 for n — —2. 
The branch that runs through the origin is associated with 
moderate density fluctuations and is the relevant one for the 
expansion (|64p in terms of cumulants at leading order in the 
quasi-linear limit. 
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FIG. 4: (Color online) The Legendre transformation (|82|l of 
the curve —t'^{Q)/2, which gives the generating function (p{y). 
The first-contact line yQ + c, of fixed slope y and height c 
decreasing from +oo, with the curve —t^{Q)/2, intersects the 
vertical axis at (0, —ip) (i.e. for c = —'p)- We show the cases 
n = and n = —2 in three dimensions. 



We display in Fig. [3] the cumulant generating function 
(p{y) that we obtain in dimension d — 3 for the two in- 
dices n = and n — ~2. This provides an iUustration 
of all the behaviors ([M|) - ([57)l . The appearance of these 
singular behaviors can also be seen from the geometri- 
cal construction of the Legendre transform ((82)) . that we 
show in Fig. [5] for these three-dimensional cases, n — 
and n = —2. For a given y, —ip is obtained as the inter- 
cept on the vertical axis of the first-contact straight line 



yG + c, of slope y, with the curve —t'^{Q)/2, decreasing 
its height c from -t-oo. Thus, the Legendre transform ([5^ 
follows the concave hull of the function —t'^{Q)/2 and it 
is regular if the latter is concave over Q < Q < oo. Note 
that this is obviously the case in the linear regime where 
T = l-g. 

For n = we have t(0) = (the curve —t'^{Q)/2 shows 
a steep up-turn at very low Q in Fig. 2]) so that for y > 
the global minimum is r = 0: the point in the range 
Qs < Q < ^ with a tangent of slope y > is only a local 
minimum and there is a local maximum in the range 
Q < Q < Qs- The local minimum corresponds to the 
regular branch in Fig. [3] and Table [TVl that runs through 
(y9(0) = 0, while the local maximum corresponds to the 
second branch in Fig. [3] and Table IIVI For n = — 2 we 
have t{Q)'^ oc Q'^/'^ at large overdensities, from Eg. ([5^ . 
Then, for < y < we again have a local minimum, 
with 1 < Q < Gs, and a local maximum, with G > Gs- 

In the quasi-linear limit, (tslt ~* such large density 
fluctuations are exponentially suppressed by a term of 
order e ''i'- , as seen in Eq. (|99|) below, so that it is 

sufficient to define the generating function by the branch 
that runs through the origin. 

4- Probability distribution V{r]r) 

Finally, from the cumulant generating function (f{y) we 
obtain through an inverse Laplace transform the proba- 
bility distribution V{r]r) in the quasi- linear limit. Using 
Eqs.linni),®, we have 

/-\-ioc J 
^[yvr-viy)]/'^!^^ ^ (98) 
-ioo 27ricr|^^ 

It is best to compute the integral ((98|l exactly, using the 
branch of (p{y) that runs through (p{y) = in case this 
function is multivalued (then it only applies to some 
range of overdensities r]r around unity). However, in 
the quasi-linear limit at fixed r/r, it is again possible to 
evaluate the distribution V{rir) from a steepest-descent 
method. The saddle-point of the exponent in ([55)1 is given 
by Vr — f'{y)- On the other hand, from ([5T|) we have 
^'(y) — G, whence G = Vr- Therefore, as expected the 
probability distribution V{rir) at point rjr is governed by 
the saddle-point described in section llV A 2l such that its 
overdensity G is equal to ry^, and we obtain from Eqs. ([5T|) . 

VM e-^^''-)^/^^^.) = (99) 

Thus, in the saddle-point approximation associated with 
the quasi-linear limit there is a precise correspondence 
between the overdensity rjr and its Laplace conjugate 
y, and the variable r expresses the Gaussian weight of 
the initial velocity (or linear density contrast) as shown 
by Ea. (P5|) . The nontrivial relation r(ryr) describes both 
the evolution ([7^ of the density of a Lagrangian region. 
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TABLE V: Asymptotic behavior of the probability distribu- 
tion Vijjr) in the quasi-linear regime (ts^^ (i.e. t — > or 
r <x), for the initial conditions of Table IH The last col- 
umn shows the range of overdensities rir where these results 
apply. If the lower threshold is not zero (i.e. —2 < n < d — 3, 
see Table HIH . it means that the spherical saddle-point forms 
shocks for lower densities. 

which only depends on the dimension d, and the effect 
([75)) . (|85p . of the change of size from g to r, which in- 
volves the initial power-spectrum index n through the 
dependence of initial power on scale. In particular, for 
the power-law power spectra P^ . using Ea. ([57|) . Eq. dMl) 
reads as 

as.^^O: lnV{vr)--:r^(v7^ . (100) 

In the large dimensional limit, d ^ oo, we have seen that 
T(r\T) = — ln(77r) (last row of Table IIV|) . whence 

In"^ ( ) 

d ^ oo, as^^ ^ : InViiir) - (101) 

OLr 

Note however that Eqs. (|100p - (|10ip only hold for densities 
T]r such that the saddle-point ((76|) has not formed shocks 
yet. As discussed in section llV A 21 and summarized in 
Table IIIIl this implies n < d — 3 and it gives a lower 
bound for ry^ if > ^2. These lower bounds are given 
by the last column in Table |V] and they will be obtained 
in section IV Bl below. 

As seen from the last expression ((99)) . the tails of the 
probability distribution V{rjr) are simply governed at 

leading order by the initial Gaussian weight e ^' ^'-i 
of the associated initial fluctuation S^q at the Lagrangian 
scale q. In fact, Eq. ([99|) could be directly obtained from 
a Lagrange multiplier method, without introducing the 
generating function ip{y). Indeed, in the rare-event limit 
we may write 

rare events : Vir]) ^ max ^-^^^--^l ■'5^. 

(102) 

That is, P(?7) is governed by the maximum of the Gaus- 
sian weight e~'^^ '-^i ■'5-l)/2 subject to the constraint 
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?7r [<^L] = V (assuming there are no degenerate maxima). 
Then, we can obtain this maximum by minimizing the 
action iS[(5l]/(t^^^ of Ea. (p5)) . where y plays the role of 
a Lagrange multiplier. This gives the saddle-point (|76p . 
and the amplitude 5Lq and the radius q are directly ob- 
tained from the constraint ?/ = J-{&Lq)i as in Eq. ([72|) . 
Then, we do not need the explicit expression of the La- 
grange multiplier y, as this is sufficient to obtain the last 
expression of the asymptotic tail ((99)) . Nevertheless, it is 
useful to introduce the generating function (/3(y), which 
makes it clear that the Lagrange multiplier y is also the 
Laplace conjugate of the nonlinear overdensity rj as in 
Eq. ((59l) . since it is also of interest by itself, as it yields 
the density cumulants through the expansion ([64]) . More- 
over, it is easier to check through the action S and the 
generating function ip{y) that the path integral (|65p is 
indeed dominated by a saddle-point. On the other hand, 
as noticed above, in the quasi-linear regime it is best to 
compute the distribution P(77r) from the integral (PS)) . 
expressed in terms of ^p{y), as the property (^(0) = au- 
tomatically ensures that the probability distribution is 
properly normalized to unity (in the case of the gravita- 
tional dynamics this has been seen to give a good match 
with numerical results [2^[3^). 

We can note that at this order the nonlinear distribu- 
tion V{r]r) could be described by a spherical-dynamics 
model, where one makes the approximation V{r]r)dr]r = 
'PLiSLq)'iSLq with rjr = J-{5Lq) and Vl is the initial 
distribution of the linear density contrast, as developed 
for instance in 55] for the collisionless gravitational dy- 
namics. Note that such a phenomenological model can 
be readily extended to non-Gaussian initial conditions. 
However, one needs the steepest-descent framework de- 
scribed in the previous sections to justify the behavior 
([99| for the rare-event tails. Moreover, in cases where 
collisions (shocks) take place, such a phenomenological 
model becomes ambiguous, while the saddle-point ap- 
proach remains valid and allows to derive exact results, 
as we shall describe in section IV Bl below. 

We show in Table |V] the asymptotic behaviors pOOp 
obtained for the initial conditions given in TablelH as well 
as the lower bound rj^ below which the saddle-point (j76p 
forms shocks. The value of 77* will be derived in sectionlVl 
below, where we take into account shocks. Note that it is 
not related to the value Qs where the cumulant generating 
function is singular, which was given in Table [TVl 

We show the distribution 7^(?/r) in Fig- El for a linear 
variance cr|^_ = 0.5 and the cases n = and n — —2 
in d = 3. For finite uslt it is better to use the inte- 
gral ([55]) . rather than the asymptotic result ([M)) . as it 
ensures that the distribution is normalized to unity and 
captures the asymmetry of the distribution with the shift 
of its peak. Note that in the case {n = —2, d = 1} the 
inverse Laplace transform (|98p of the quasi-linear gener- 
ating function given in the fourth row of Table IIVI gives 
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d = 3 




FIG. 5: (Color online) The probability distribution Vijjr) in 
the quasi-linear limit from Ea. (|98p . We show the cases n = 
and n = —2 in three dimensions, for a linear variance a^^ = 
0.5. 

the explicit expression 

n = -2,rf = 1 : VM = e-(^-7fe)'/(^<.-\ 

(103) 

Again, as seen in [l^l, it happens that in this case the 
result (|103|) is actually exact. In generic cases, deviations 
from the quasi-linear limiting distribution should appear 
at small scales and late times. 

From the geometrical construction described in Fig. 3] 
we can see that singularities for the function ^p{y) occur 
at inflexion points of the curve — r^(tJ)/2. In particular, 
in agreement with the analysis of Eqs. ([92l) - (f94|) . at large 
densities a singularity appears as soon as /Q for 
Q oo, which implies that —t'^{Q)/2 is no longer con- 
cave at large Q. From Ea. ([M)) this simply corresponds to 
a sub-exponential large-density tail of the form ~ e"''-- 
with a < 1. Then, it is clear that the integral ([55]) is 
divergent for y < so that the exact cumulant generat- 
ing function has a branch cut on the negative real axis, 
y < 0, even though the cumulants of all orders may be 
finite. In this case, the singularity at < and the 
two branches observed on yg < y < are related to this 
branch cut and to this sub-exponential large-density tail. 

For the low-density tail, this construction shows that 
when t(Q) and t'{Q) remain finite for Q —^ 0, a, singu- 
larity appears on the positive real axis, ys > 0, or ip{y) 
is restricted to a finite range y < ys- Thus, singularities 
on the positive real axis are generically associated with 
distributions that do not vanish in the limit rjr —^ 0. 

In any case, in the quasi-linear limit the large- 
fluctuation regime associated with the second branch of 
f{y) is irrelevant, as it is exponentially suppressed by a 

factor of the form e ''^''lt^ see Eq. ([M)) . Therefore, we 
only plot in Fig. [5] the range associated with the quasi- 
linear branch of fiy), that is, r]r < 2\/2 for n — —2, 



see Table HVl For n = 0, this would require r]r > 1/8, 
but as shown in Table IVl and explained in section fV B 71 
below, shocks appear before this threshold, as soon as 
ijr < (2/3)'^, below which Eq. (|100p is no longer valid. 
Therefore, in the case {n — 0,d — 3}, we only plot 
the result ([98]) . obtained from the quasi- linear generating 
function, above this lower-density bound, rjr > (2/3)'^. In 
the limit cr^^^ the weight associated with such low- 

density regions decays exponentially as e ' ^i-^- (disre- 
garding the numerical factor), but at a^^ = 0.5 this 
region is already non-negligible. Note that the upturn at 
low density shown in Fig. [5] is not necessarily a signature 
of the breakdown of the quasi-linear limit in this domain 
at ag^ =0.5. Indeed, as discussed in section IV B 7[ the 
distribution V{rir) does not decay exponentially at low 
density, see also Table IIXI a nd it may even grow as a 
power law. As shown in [30|, this is for instance what 
happens for the case {n = 0,c? = 1}, where at low den- 
sity the probability distribution shows an inverse square 
root tail, P(?/r) fx l/y/fj^. 

B. Velocity divergence (i.e. spherical velocity 
increment) 

We now consider the probability distribution, V{Qr), 
of the mean velocity divergence 8^, defined in Eq. ([2T|) . 
that is also the velocity increment over distance 2r av- 
eraged over all directions (up to a normalization factor). 
Following the method described in section IIV Al for the 
spherical overdensity 77^, we can also obtain the quasi- 
linear limit of 7^(8^) by a steepest-descent approach. 
Thus, we again define the moment and cumulant gener- 
ating functions '^{y) and (p{y) as in Eas. (|59|) - ([64)) . which 
can be expressed as a path integral (|65p with an action 
as in (j66p . where the nonlinear functional Tyrii^L] 
is replaced by 8r[i5L]. The action still being spherically 
symmetric, we can also look for a spherical saddle-point. 
If particles on the sphere S have not been shocked yet 
they have kept their initial velocity Uq, so that 8^ only 
depends on the initial radial velocity Uq^ of the particles 
at the Lagrangian radius q that have moved to radius r 
at time t. Then, by the same reasoning as for the over- 
density, we obtain as expected the same profile ([75)1 for 
the saddle-point, but a different function T{5Lq)- Indeed, 
since the radii q and r are related hy r — q -\- 1 uoq we 
have from Eq.(pT|) 

Or = ~^tUr{t) = ~^tUaq = d(^'^ -ij. (104) 

Note that Eq. (|104p shows that for spherical dynamics the 
quantity 8^ obeys 

8^ > -d. (105) 

This lower bound corresponds to very fast expansion, so 
that the particles observed at radius r come from the cen- 
ter q = 0. Note that this requires the initial velocity field 
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to be singular at the origin, whence n > —1, in agree- 
ment with the exponential falloff (|129p . (|134p . obtained 
at low densities for n < — 1, as shown in section |V] below. 
In such cases, where rarefaction regions can appear (i.e. 
truly empty regions), it may happen that the cell of ra- 
dius r is within a larger empty domain, so that there are 
no particles on the sphere S. However, in the quasi-linear 
limit, where we consider small values of |6r|, we do not 
consider this case. Next, from Eq. pO]) the linear density 
contrast is given by 



whence 

= J^{SLq) with J^iSLq) = -d 



1 



d 



(106) 



(107) 



Again, since the function T simply describes the spherical 
Burgers dynamics it only depends on the dimension d 
and not on the initial power spectrum. Then, we obtain 
the expression (j78p with this new function JF, and we 
can define the associated variables r and Q as in (|79p . so 
that the cumulant generating function (p(y) is given by 
the relations ((8T|) and ([82]) . The spherical saddle-point 
being identical to the one obtained in section llV A21 the 
profile (|83p still applies for power-law power spectra, as 
well as Figs. [H El In particular, if n > d — 3 shocks 
appear as soon as t 7^ so that this saddle-point is no 
longer exact in such cases. As in section IIV Al we shall 
not consider the modifications that appear in such cases, 
as they should remain small in the quasi-linear regime 
(with d— 3 < n < d—2), and we focus on cases such that 
n < d — 3. 

Close to the origin, since by symmetry we have (u) — 0, 
whence (6r) = 0, we always have (compare with ([84]) ) 



0, y^O: Q T, y 



and 



y 

f ;r, 



(108) 

and keeping only these low-order terms gives back the lin- 
ear Gaussian of variance cr|^ . For the power-law initial 
power spectra ([6]) we obtain from Eqs. ([79| . (|104p . 



which leads to 



~d 



T{g) = d 1 + ^ 



(109) 



(110) 



and to the parametric representation of ^{y)^ 



-dy-\ 



n+l 



2{n+l) 1 + ^ 
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(111) 
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FIG. 6: (Color online) The cumulant generating function ifijj) 
for the velocity divergence Qr, in dimension d = 3 for the 
power indices n = (solid line) and n = —2 (dashed line), 
from rows 2 and 4 of Table IVII In both cases there is a 
singularity on the real axis, with ys — 1/2 for n = and 
ys — —3/2 for n — —2, but only one branch for n = —2. 
The branch that runs through the origin is associated with 
moderate velocity fluctuations and is the relevant one for the 
expansion in terms of cumulants in the quasi- linear limit. 



2\ d 



(n+1) 1 + 



n+l 



^\ n+2 / ^\ n+3 

-2(n + 2)(l + ^j +(n + 3)h+^ 



(112) 



We can see that, contrary to the Eas. ((55)) - ([M)) associ- 
ated with the density contrast, the dependence on d of 
the system (|llip - ()112p simplifies as d~'^(p{dy) no longer 
depends on the dimension d. This implies for the proba- 
bility distribution the scaling 



Or 

~d 



d2 



(113) 



where we noted 75^(6^; erf ) the quasi-linear probability 
density of in dimension d when the linear variance is 
CT^^^. Thus, in the quasi-linear limit the change of di- 
mension is fully absorbed by a rescaling of the velocity 
divergence Qr and of the linear variance crf^^. There- 
fore, contrary to the case of the overdensity studied in 
section ITV Al the properties of (p{y) and V{Qr), such as 
the presence of singularities and sub-exponential tails, 
only depend on n and not on the dimension d. 

Expanding near the origin, we obtain for the third and 
fourth-order cumulants in the quasi-linear limit: 



0: & 
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= -3- 
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3(n+ l)(7?i + 9) 

d2 



(114) 



(115) 
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— oo < y < d 


— d < CJ < oo 








-2 




-f < y < oo 


— d < C/ < oo 








d = oo 


2 


— oo < y < oo 


— oo < Q < oo 









TABLE VI: The cumulant generating function (p{y) of the velocity divergence O^, in the quasi-linear limit (t|^^ — > 0, for integer 
values of n and arbitrary dimension d, where explicit solutions of the system piH) - (|112|) can be obtained. The columns 2-4 
show the quasi-linear branch, and the associated range of {y,G} that contains the point {0,0}. The columns 5-7 show the 
second branch, associated with very rare events, that appears in some cases (here only for n = 0). The last row shows the 
infinite-dimensional limit, d oo, which no longer depends on n and corresponds to the Gaussian. 



As discussed below (|2T|) this also gives the cumulants 
associated with the third and fourth-order spherical ve- 
locity structure functions in the quasi- linear limit. As 
for Eqs. ((90|) - ((9T|) . these quantities are not the standard 
skewness and kurtosis, and the powers in the denomina- 
tors are such that they have a finite non-zero quasi-hnear 
limit. 

For integer values of n we can also derive explicit solu- 
tions to Eas. pTT|) - pT2)) . which we show in Table IVll As 
for the density studied in section HV A 3[ the quasi-linear 
generating function ip{y) can show two branches when 
the function y{G) is not monotonic over G €] — d, -|-oo[ 
(the variable G now covers the range ] — d, +oo [, as seen 
from Eq. (|105p ). This occurs for n = 0, shown in the sec- 
ond row in Table IVTl while for n = — 1 and n — —2 there 
is only one branch. 

In dimension d = 1 we have from Eas. (|104p . ((72|) . 



in Table EI 



Vr 



(116) 



so that the distributions 'Pi{r]r) and 7^1(9^) are identical 
up to a shift of unity. Then, we can check that the results 
shown in Tables IIVI and IVII for the case {n — —2, d — 
1} are consistent. In particular, the generating function 
given in the fourth row of Table IVll yields the probability 
distribution 



V{0, 



X exp 




(117) 



We can check that this agrees with relation (|116|) and 
Eq. (|103p for d = 1. Moreover, in dimension d = 1 this 
result (|117p again happens to be exact [Tsj . 

For n = — 1 we simply obtain the Gaussian (third row 



(118) 



Thus, the effects of the nonlinear evolution (|107p and of 
the change of scale q ^ r (encoded in the change from 
J-' to G) compensate in such a way that at leading or- 
der in the quasi- linear limit the cumulants (6r)c vanish, 



whence (8J!)c ^ cr, 



2(p-l) 



for asr 



and p > 3, in 



agreement with Eqs. (lll4p - (|115p . However, note that the 
distribution (|118p differs from the linear Gaussian in the 
sense that is restricted to 8^ > —d, from Eq. (|105p . 
Of course, the weight associated with this lower bound 
becomes exponentially small in the quasi-linear limit, so 
that it cannot be seen in the leading-order value of the 
cumulants (0^)c, whence in the quasi-linear generating 
function ip{y). 

On the other hand, in the limit of large dimension, 
d —f oo, we obtain tj(r) = — r which gives back the linear 
Gaussian of Eq. ljllSp . again in agreement with Eqs. (|114p - 
(jllSp . However, contrary to the distribution (jllSp asso- 
ciated with n = —1 at finite d, in the limit d — > oo the 
Gaussian extends down to — oo, since the lower bound 
pOSp is repelled to — oo. Note that this was not the case 
for the overdensity rjr, where the probability distribution 
did not tend to the Gaussian but to a lognormal dis- 
tribution for d — > oo, see the last row in Table IIVI and 
Eqs.(inni)-(EI1)- 

We show in Fig.[5]the quasi- linear cumulant generating 
function ip{y), obtained in three dimensions for n — and 
n = — 2. Although there is a singularity on the real axis 
for both cases, there is only one branch for n — —2. 

Next, the quasi-linear probability distribution V{Qr) 
is obtained from the cumulant generating function f{y) 
by an inverse Laplace transform, as in Eq. (|98p . In the 
quasi-linear limit, as for the overdensity rjr, it obeys the 
asymptotic behavior (|99p. Then, using Eq. (|110p we ob- 



20 



n 


d 




Br 





3 




2 


e, > -1 


-1 


2 




> -1 


-2 


1 


[(i+e.)"''" - (1+e,.)''" 


2 


> -1 


-2 


3 




r 


Or > -3 


n 


oo 







TABLE VII: Asymptotic behavior of the probability distribu- 
tion ViQr) in the quasi-linear regime, for a few integer values 
of n and d, using the normalization of Table |l] for the initial 
conditions. The last column shows the range of spherical ve- 
locity increment Or where these results apply. If the lower 
threshold is not equal to —d (i.e. —2 < n < d — 3), it means 
that the spherical saddle-point forms shocks for lower Or. 
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FIG. 7: (Color online) The probability distribution V{Or) of 
the velocity divergence Or (spherical velocity increment) in 
the quasi-linear limit. We show the cases n = and n = —2 
in three dimensions, for a linear variance (t| =0.5. 
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- 1 
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and in the large dimensional limit, where t — —Or 
d ^ oo, (75^^ — > : lnV{Qr) 



(119) 



Ol 



(120) 



which only hold for n < d—3 and above a low-O^ thresh- 
old if n > -2, see Table Uni Note that Eq. (fTT9l) may be 
directly obtained from the asymptotic tail (jlOOp derived 
for the density distribution, V{r]r), by substituting 



d(^y''-l), r^^^^l + ^y. (121) 



This relation follows from Egs. fT^ . p04p . that express 
both rjr and in terms of the Lagrangian radius q of 
the saddle-point, so that there is a unique correspondence 
between rj^ and Or for these spherical saddle-points. 

We show in Table IVIII the asymptotic behaviors pi9p 
obtained for the initial conditions given in TablelH as well 
as the lower bound G* below which the saddle-point ([76]) 
forms shocks, which will be derived in section |V] below 
where we take into account shocks. Again, the value (and 
the existence) of 9* is not related to the value (and the 
existence) of Qg where the cumulant generating function 
is singular, which was given in Table IVTl 

We show in Fig. [7] our results for 7'(8r) for the cases 
n — and n = —2 in three dimensions, for a linear 
variance a^^^ = 0.5, as in Fig. [5] Again, we only plot 
the distribution over the range associated with the quasi- 
linear branch of (p{y). For n = —2 this actually covers the 
whole range Qr > —d, but for n = this corresponds to 



Qr > —d/3, see Table lVTl Note that in the latter case this 
lower bound happens to coincide with the lower bound 
0* where shocks appear (for c? = 3). The comparison 
between Figs. [7] and [5] shows that the distribution of the 
velocity divergence 0^ remains closer to the Gaussian 
than the distribution of the overdensity rjr- This can also 
be seen from the fact that the singularities tJs are farther 
from the origin y ^ (compare Tables I VII and IIVI and 
Figs, ini and [3]). On the other hand, in agreement with 
(|114p . we can see that the skewness has opposite signs 
for n = and n = —2, as the peak of the distribution 
shifts to either side of 8r = 0, while the skewness of the 
density was always positive for n < d — 2, see Eg. ([90)1 . 
which covers the range where cr^^^ is well defined. 



ASYMPTOTIC TAILS 



The results obtained in section Ivl applied to the quasi- 
linear limit, as^r ~^ 0? for the case n < d — 3, so that 
as^r is well defined and shocks only appear after a finite 
time (if d — 3 < n < d — 2 shocks appear as soon as 
t but CTaj^^ is still well defined and these results should 
provide a reasonable approximation). We now consider 
the limit of rare events, that is very large density and 
velocity fluctuations at flxed linear variance a^^ , or at 
fixed (7^^^ if ag^^ is divergent. Thus we study the tails of 
the probability distributions 'P{rjr) and 7^(0^), for any 
value of cfs^^ or a^p^^. As in section [IVI we can use a 
steepest-descent approach and look for the minimum of 
the action 5, defined as in This will give the tails 

of the probability distribution through Eq. ([M|) . 
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A. Saddle-point without shocks 

1. Rare density fluctuations 

For n < d — 3 we can use the same action S[Sl] as 
in (|66p and we obtain the same saddle-point defined by 
Eas. ((76)l and ((82|) provided no shocks have formed. As 
discussed in section IIV A21 this constraint is satisfied 
for overdensities if n < d — 3 and for underdensities if 
—3 < n < —2. Therefore, in such cases Ea. ([M)) re- 
mains vahd, where T{r]r) is still given by Eq. ((87|) . pro- 
vided the saddle-point approximation is justified. Thus, 
we must show that as we consider very large density fluc- 
tuations the Laplace transform ^(y) and the distribution 
V{rir) are dominated by an increasingly narrow region 
around this saddle-point. We no longer have a fixed ac- 
tion S[Sl] multiplied by a prefactor that goes to infinity, 
as was the case in the quasi-linear regime for (I65p . with 



oo. Therefore, the analysis is more compli- 



cated as we should study the Hessian of the action at 
the saddle-point, which requires second-order perturba- 
tion theory around this saddle-point, taking into account 
angular degrees of freedom. Here we shall simply show 
that the steepest-descent approximation is well justified 
for ^'(y) (i.e. <^(y)) with respect to the family of ini- 
tial states ((76)) - (|77l) . parameterized by the overdensity Q 
within radius r. 

In this subspace, the action S[Sl] is an ordinary func- 
tion SiG), given by Eq.® as SiG) ^ yG + t{G)^/2. 
Thus, at the saddle-point the first and second derivative 
read as 



dG 



0, 



d^ 
dG^ 



dG'^ 



r{Gf 



(122) 



Let us first consider the case of large underdensities, G 
0, with -3 < n < -2. Then, from Eq.lUZl) we have 



: 



d^^ (n + l)(n-|-l-d) 
W 2 



g— 



(123) 



Then, over this subspace, disregarding prefactors asso- 
ciated with changes of variables, we write the analog of 
([55)1 as 



dGe 



-S{Q)/o 



(124) 



and expanding the action around the saddle-point Gc^ 
defined by the condition dS/dG = 0, we find that only 
values sufficiently close to Gc contribute to the integral, 
with 



\G-Gc 



1 



(125) 



whence (since —3 < n < —2) 

g^^O: ^ g-(n+i)/(2d) ^ 0^ (^26) 



Since S{g) behaves as a power law, Eq. (|126p implies that 
higher-order terms beyond the Gaussian around f/c are 
subdominant. Thus, at large positive ?/, which corre- 
sponds to C/c — > 0, see (inH), the Laplace transform tf{y) 
is dominated by its saddle-point (if we only consider the 
subspace described by the profile ([76]) ). Then, assuming 
that this remains true when we take into account other 
degrees of freedom, and that there is no deeper minimum 
associated with strong deviations from spherical symme- 
try, we recover at large y the behavior (j95p obtained in 
section IIV A 31 where we studied the quasi-linear regime. 
Next, the distribution V{r]r) is still given by the inverse 
Laplace transform (|98p . and it only remains to show that 
this integral is again dominated by its saddle-point, in the 
limit rjr — > 0. As noticed in section IIV A 41 the saddle- 
point yc of the exponential satisfies rjr = Lp'{yc), whence 
rjr — t/c where is related to yc through the Legendre 
transformation (|8ip . Thus, in the limit 77^ — > we ob- 
tain j/c — ^ +00 and from the previous results the behavior 
([55)1 . Next, expanding the exponential around y^ we find 
that contributions to Virjr) come from the range 



\y-yc\ 



1 



yc 



' 2(n + l-d) 



where we used the large- ?/ behavior ([55]) . This yields 



(127) 



yc 



00 



\y-yc\ 
yc 



yc 



0, 



(128) 



which again ensures that the integral ([98)1 is dominated 
by the Gaussian integration around its saddle-point in 
the limit rjr 0. Therefore, the steepest-descent ap- 
proximation is legitimate at very low densities, whatever 
the value of as^^ , and we recover the rare-event tail (|99p . 
As in (|100p . this yields the low-density tail 



3 < n < -2, 7]r ^0: VM 



(129) 

We can apply the same procedure for large overden- 
sities, rjr +00. Here, a subtlety arises from the fact 
that for n < d — 3 large densities are associated with the 
second branch |[M|) of (p{y) discussed in section IIV A 31 
As recalled in section IIV A 41 this simply corresponds 
to cases where the high-density tail of the distribution 
V{rir) shows a sub-exponential decay, such as e"''" with 
a < 1. From the definition ([55)1 we can see at once that 
in such cases the Laplace transform ^'(j/), whence f{y), 
diverge for negative y: they generically have a branch cut 
along the negative real axis. Then, the saddle-point in 
the path integral (j65p becomes a local maximum but one 
can still apply the steepest-descent method, using an ap- 
propriate deformation of the integration contours in the 
complex plane. This is discussed in details in section 3.6 
and appendices A and B of ^ hence we do not give 
further comments on this point here, see also the chapter 
on instanton contributions in [4^ . Then, as for under- 
densities, we can check that for large overdensities the 
contributions to f{y) come from an increasingly narrow 
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range of Q, when we restrict to the subspace parameter- 
ized by the profile ((75|) . as in (|124p . Thus, using again 
Eq. ([87|) , we now find that the range of overdensities that 
contribute to (p{y) behaves as 



IG-Gc 



g-(n+3)/(2d) _^ 0^ (^30) 



SO that (p{y) is still governed by the saddle-point (as- 
suming there are no larger contributions associated with 
strong deviations from spherical symmetry), and we re- 
cover the behavior ([94|) . Then, the range of y that con- 
tributes to the large-density tail is 

Vc^O-: k-JlA ^ \y^\2U^ ^ 0. (131) 

\yc\ 

Therefore, the steepest-descent approximation is again 
legitimate at very high densities, whatever the value of 
as^^, and we recover the rare-event tail ((99)) . As in (|100p . 
this yields the high-density tail 



case particles that would reach the center would stick 
there. Thus, the infinitesimal viscosity regularizes the 
dynamics and makes the sensitivity to non-spherical per- 
turbations finite. This well-behaved dependence on the 
initial conditions clearly appears through the Hopf-Cole 
solution recalled in section Hi Dl 

If we are only interested in the exponents that appear 
in the expressions (I129P , (|132p , disregarding the numeri- 
cal factors and using Eg. ([29]) they read at leading order 
as 



n < d — 3, Tjr ^ 
-3 < n < —2, rjr 



: \nV{r]r) oc - r''+^ T]f^ jt^ {IM) 



We give in Table lylTTl the explicit expressions of the tails 
(|129p . (|132p . for the initial conditions normalized as in 
Table [H and we mark as "shock" the cases where the 
saddle-point discussed above gives rise to shocks. This 
agrees with Table|V]where we only keep the leading term 
for r/r CO or rjr 0. 



n<d-3, rir^oo: VM - e~'''''-'^/^^'"'i.-\ (132) 

We can check that in the case {n = —2, d — 1} both 
asymptotics (|129p . (|132p . agree with the exact distribu- 
tion (fTU5)) . 

In the large dimensional limit, d — > c», we still obtain 
Eg. ljlOip . which applies to both rare overdensities and 
underdensities, but only if n < —2 in the latter case. 

In fact, the saddle-point approximation is valid as long 
as we consider the limit of rare events, which corresponds 



either to the quasi-linear limit, (JSl,- ~* (i-^- 



oo) at fixed density rjr and velocity increment 0^, or 
to the limit of extreme densities, 77^ or 77^ — s- oo, 
and extreme velocities, Qr —d or — > 00, at fixed 
time and scale (i.e. at finite cri^^). Of course, the range 
of density fluctuations to which these results apply is a 
priori repelled to increasingly large fluctuations, rjr — > 
00 or 0, as cfsl,- grows (since they correspond to 

rare events, V 1). In such regimes, integrals such as 
(I65p are governed at leading order by the minimum of 
the action 5, which shows a steep dependence on the 
initial conditions. This legitimates the steepest-descent 
approximation in these cases. 

We can note that an exception to this behavior occurs 
when the action happens to be singular at its minimum, 
that is the initial conditions where S is close to its min- 
imum form a subspace of measure zero. For ordinary 
integrals such as p24p . this corresponds to cases where 
the function S{G) is discontinuous at the point Gc^ with 
S{Gc) being strictly smaller than both left and right lim- 
its, S{Gc) s-iid 'S{G^)- This possibility actually appears 
for the case of the collisionless gravitational dynamics 
(sgI , where at large positive densities a strong radial orbit 
instability appears (associated with the extreme sensitiv- 
ity of the trajectories, that actually diverges, as particles 
move through the center of the object). This problem 
does not appear in the Burgers dynamics, since in any 



2. Rare velocity fluctuations 

The discussion oiV{rir) in the previous section directly 
extends to the tails of the distribution V{Qr) of the mean 
divergence Qr (spherical velocity increment), deflned in 
(PT|) . Thus, for the same cases as in section FV A II where 
the saddle-point does not form shocks, the rare-event 
tails are still given at leading order by the flrst exponen- 
tial in Eq. ([g5|) . where T(9r) is given by Eg. piOp . whence 
by Eq.dnil). This yields 



< d - 3, 



d^ fQr 



2aj 

OL 



n+3 



OC 



i2 



(135) 



and 



3 < n < -2, e,. 
InP(e^) 



> -d : 
-d^ 

OLr 

^n+3 



n+1 



OC 



i2 



1 



d 



n+1 



(136) 
(137) 



In the large dimension limit we still have Eq. p20l) . which 
holds for Qr 00 for any n, but for 0,. — s- —00 only if 
n < — 2. We also give these results in Table Ivnll 

Note that the appearance of a second branch for (f{y) 
is not simultaneous for rjr and Qr, as noticed in sec- 
tion IIVBI and as can be seen from the comparison of 
Tables IIVI and IVII Since for spherical initial conditions 
rjr and 8^ are related by (|12ip . it is clear that the sub- 
tlety associated with the second branch of ip{y), that is, 
the branch cut of the exact Laplace transform, is only 
a mathematical difficulty due to a sub-exponential tail 
but has no physical effect. Thus, in cases where only 
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n 


d 


InVirjr) for rjr ^ oo 


ln'P{rjr) for 77^ ^ 


InP(er) for Qr ^ 00 


lnV{Qr) for e,. -d 





3 




shock 


-r^e?/(54t2) 


shock 


-1 


2 




shock 


-r2e?/(8t") 


shock 


-2 


1 






-re,./(2t2) 


-r/\2t\l + Qr)] 


-2 


3 




-5r/(8i'r?y') 


-5rer/(24f^) 


-15r/[8t^(3 + er)] 


n 


oo 


^-^72 In (rir) 


(241^ In^ (,,,,) ifn< 2 


, 1 + 3 

2fi '^r 


, 1+3 

(2r^) 2 „2 -r ^0 
— -^^ 2*2 11 < — 2 



TABLE VIII: The tails of the distributions Virjr) and P(0r) for a few integer values of n and d where the relevant saddle-point 
is regular (i.e. does not induce shocks), from Eqs. p29|) . (|132p . The initial conditions are those given in Table |I] (i.e. with the 
same normalization). Cases marked as "shock" correspond to saddle-points that give rise to shocks, so that the normalization 
factors in Eqs.([T29}, (|132|) . and (|135p - (|136p . are no longer valid. They are analyzed in section FV Bl (this includes the case 
{n = 0, d = 1} for both large and small densities). For d — > 00 the expressions given at low densities and velocity divergences 
only hold for n < —2 {n > —2 gives rise to shocks). 



the density quasi-linear generating function ip{y) shows 
a second branch, we can first compute the high-8r tail 
of 7^(0^), which only involves local minima of the ac- 
tion and does not require a deformation of integration 
contours, and next use the relation (|12ip to derive the 
high-density tail of V{rjr). More generally, as noticed in 
[39|, in order to avoid the complications associated with 
sub-exponential tails, we can simply compute the rare- 
event tail of the quantity X = rj^ , with a small enough 
f3 so that ViX) shows a super-exponential decay, and 
next derive 'P{rjr) from V{X) through a simple change 
of variable. On the other hand, we could directly ob- 
tain Eqs.dngi-dnTl) from Eqs.dnS]), by using the 
relation (fT2T|) . 

B. Saddle-point with shocks 

1. Paraboloid construction and action S[il}o\ 

For cases where the constraints in (I129p , p32p , are not 
satisfied, that is when shocks cannot be ignored, we can 
no longer rely on the regular saddle-point of the previous 
section lV Al However, we can use the exact solution 
and the geometrical construction ([55]) , to study solutions 
of the equations of motion that contain shocks. In partic- 
ular, following the approach used in the previous sections, 
we can look for saddle-points of an appropriate action, 
that include shocks, to obtain the tails of the distribu- 
tion V{rir). However, to use ([36]) it is more convenient to 
work with the velocity potential tpo , rather than with its 
Laplacian, 60 = SL/t, that we used in section HVl Thus, 
we now write the Laplace transform (|59p as 

*(y) = [detC^y^^ Jv^Po e-y^-i^^^-i^^-^^o-'^o^ 

(138) 

where 77r['0o] is now the functional that affects to the 
initial condition ipQ the nonlinear overdensity rjr, built at 



time t within the cell of radius r centered on the origin, 
X = 0; and C^o(xi,X2) is the two-point correlation of 
the initial potential. Although this is no longer essential 
(since we do not consider here the limit a^„^ 0) we 
rescale the transform ^'(y) in a fashion similar to Eg. ([55)1 . 

■^(y) ^ e-^^y<o.'>/<o., (139) 

so as to obtain an action iS[7/'o] that is similar to Ea. (|66p . 

2 

s[%] = yvAH + ^ i'o-C^'-i'o, (140) 

where ct^q^ is the variance of the initial radial potential 
tpor at radius r, defined in Eq. (p4|) . 

As in section IIVI we can look for spherical saddle- 
points. Then, we can look for the minimum of the action 
5[?/'o] within the subspace of spherically symmetric ini- 
tial conditions, where "00(1) = ^09 with q ^ |q|, and the 
restriction of the action to this subspace reads as 

2 

5[^o,'] = yVrM + ^ V^,-C^J„.^o,^, (141) 

where C^g^{q[, is the radial covariance introduced in 
Eq.®. 

2. Regular saddle-point without shocks 

Let us first check that in the case where there is no 
shock, we recover from (|14ip the results obtained in sec- 
tion |TV1 As in Eq.dni), a saddle-point of the action (|14ip 
is characterized by 

V'o,' = ^ rdq"C^,M',<l"):^- (142) 

Following the discussion below Eg. ([75)1 . since 7/^ only 
depends on the initial velocity at the Lagrangian ra- 
dius q, and UQq — —dipoq/dq, the functional differential 
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Dijr/'Dipoq'i is zero for q" ^ q. However, it is no longer 
a Dirac, S^iq" — q), but its first derivative, 5'^{q" — q). 
Indeed, from the geometrical construction ((38|l . the La- 
grangian radius q is obtained as the first-contact point of 
the paraboloid 'Px,c(q) with the initial potential "iAolq)- 
Using the spherical symmetry this corresponds to the 
first-contact point of the parabola Vr.c{q) with the curve 
■009, which is characterized by the two equations (for con- 
tact and tangent slopes) 



(g - rf 
2t 



+ c, 



Oq 



(143) 



Then, as we change the initial potential by an infinites- 
imal perturbation Aipoq', both c and q are modified by 
amounts Ac and Aq, and the second Ea. (|143p gives 



= ( ^ - 



Oq 



(144) 



Then, since ry^ — (q/r)'^ as in Eg. ([72]) . we have Ar/^ « 
Aq (X AipQq, whence 



•Oq" 



^6'o{q"-q). 



(145) 



Thus, in agreement with the previous discussion, the 
derivative (jl45p vanishes for q" ^ g, but it is now the 
first derivative of the Dirac distribution. Substituting 
into Eq. (|142p gives the linear profiles of the saddle-point 
as 



d 

V'o?' oc —C^,^,Xq\q), 



(146) 



dq'dq 



Ci,oAq, q) oc Cuo^iq, q), (147) 



UOq' OC 



where we used Eq.([32|). Then, the comparison with 
Eg. lfZT)) shows that we have obtained the same spheri- 
cal saddle-point as in section lTV A 21 This means that we 
recover the results of section IIVI in the quasi-linear limit 
(here cr^or ~^ 0), and of section FV Al in the appropriate 
rare-event limits (rjr — > or 77^ — ^ 00). In particular, the 
tail of the probability distribution reads at leading order 
as 



(148) 



where the exponent is evaluated at the saddle-point, as 
in Eq. (|99|) to which it is equivalent. 



3. Taking shocks into account 



The advantage of the formulation (|14ip in terms of 
the potential ipQ is that we can now handle cases where 
shocks must be taken into account. Note that this applies 
in particular to the cases n > d — 2 where the variance 
ag of the linear density contrast is divergent. To this 



order, we generalize the previous configuration, with a 
unique first-contact point between Vr,c{q) and ^oq, to 
states where the initial potential follows the parabola 
over a finite range [g_,(7+], and remains below it else- 
where. This corresponds in particular to a shock at ra- 
dius r that contains all the matter that was initially lo- 
cated within Lagrangian radii q- and g+ (all this matter 
merging at position r at time t). Then, the nonlinear 
overdensity 77,- is not modified by infinitesimal perturba- 
tions At/jqij over q ^ q+\ (since they do not affect the 
first-contact parabola) and Ea. (|142p implies 



g > : Voq 



9+ 



dg'Cv,o,,(g, q)f{q), 



(149) 



with some kernel f{q') to be determined (in the case of a 
unique contact point, i.e. no shock, we have seen above in 
Eq. (|146p that we have f{q') oc (5^(g' — g+) and q- = q+). 
Let us note C^q^ the restriction of the kernel C^g^{qi, 52) 
to the interval [q-,q-f-]. Then, since V'Og = 'Pr,c{q) over 
the range Ea. (|149p implies over this interval: 

q- <q< q+ ■ Vr.c = C^^^^.f, whence / = C^^ .'Pr,c 

(150) 

Moreover, substituting Eq. (ll49p into the action (|14ip 
yields 



S = yrir 



war 



dgi dq2 / (<7i ) C^or (91 : 92 ) / (92 ) (151) 



9+ 



--yVr + ^J d(?idg2Pr,c(<Zl)C0„^(<Zl,g2)Pr,c(g2)(152) 



a-, ri+ 



9- 



dqf{q)VrAq)- 



(153) 



^ -1 



Next, since rjr and (T^t^only depend on q_ and min- 
imizing the action (|152p with respect to the height c of 
the parabola Vr,c gives 



dS 



1+ 



— = / AqiAq2 C^^^iqi,q2)VrAq2) - 0, (154) 



whence, using Eq. (|150p . 



9+ 



dg/(q)=0. 



(155) 



Thus, in order to minimize the action iS['0o] over the 
spherically symmetric initial conditions that show a 
shock at radius r we proceed in two steps. First, the min- 
imum of S over the class of profiles ipQ that follow their 
first-contact parabola over a finite range [q^,q^], to be 
determined afterwards, is obtained by solving Eqs. (|150p 
and (jl55p . This gives the kernel / and the parabola 
height c as a function of the parameters g_ and (7+. 
Second, substituting into the action (|153p we minimize 
S over q^ and q^. This provides the minimum of the 
action p4ip over all spherical states such that the Eu- 
lerian radius r maps to a continuous Lagrangian range 
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q+], provided the saddle-point obtained in this fash- 
ion remains strictly below the parabola Vr^c outside of 
the interval which we must check afterwards 

since we have not imposed the constraint ipQ < Vr,c hi 
the previous derivation. A priori it could happen that 
the minimum of the action is reached for initial config- 
urations that touch the first-contact parabola over dis- 
joint regions. Then, this would be seen by noticing that 
the minimum obtained through the previous procedure 
touches or crosses the parabola P^.c somewhere outside 
of the range In such a case, one would need to 

generalize the approach described above to initial states 
that follow their first-contact parabola over k several dis- 
joint intervals [q^^ ,q^_^], i — l,..,k. Starting with the 
case k = 1 discussed above, one could add a new contact 
interval in a series of steps, until the minimizing pro- 
file remains below the first-contact parabola everywhere 
outside of the k contact intervals. Note that this method 
also includes the case where some intervals are reduced 
to a point, which corresponds to the limit q^_^^ — q^l^ — > 0. 
Thus, this covers the case where the parabola would only 
have two (or a few) isolated first-contact points. 

Finally, we must specify the relation between the over- 
density rjr and the Lagrangian coordinates q- and g+. 
We must separate the cases of large overdensities and 
underdensities as 



?7r > 1 : Vr 



r]r <l 



n. 



(156) 

Indeed, since we have a shock at radius r, with a finite 
mass TTT.'^hock ^ _ ^Yie density within radius r and 
the enclosed mass are ambiguous. It is actually discontin- 
uous at r, going from m_ to m+ with m+ — m_ = m^^°'^^. 
This also means that the minimum discussed above is 
unstable, in the sense that an infinitesimal perturbation 
will move this mass inward or outward, so that m{< r) 
goes to TO_ or m+. Then, if we consider rare and large 
overdensities, the probability 7^(?/r) will be governed by 
the initial conditions close to the previous minimum such 
that m(< r) ~ m+, which leads to the first Eq. (|156p . 
Similarly, for extreme underdensities we obtain the sec- 
ond Eq. (|156p . Note that this also shows that the action S 
is not regular at the minimum obtained above and going 
beyond the leading term given by the exponential as in 
Eq. (|148|) would require a careful analysis. However, this 
discontinuity is not of the same kind as the one encoun- 
tered for coUisionless gravitational collapse, associated 
with radial orbit instability and recalled in section IV A II 
above, since by using (|156p we simply consider the left or 
right limit of S (with respect to any degree of freedom) 
and not an isolated point of zero measure. 



4- Computation of the saddle-point 

In practice, it can be difficult to solve Eg. pSOp for the 
kernel /, and we did not obtain a general solution. How- 
ever, for the power-law power spectra ^ with low integer 



values of n and d, where the radial potential correlation 
C^fto,, (<7i, (72) takes the simple forms given in Table IH it 
is possible to derive explicit expressions for /(g) from 
Eq. (|150p . It is convenient to first write f{q) as a deriva- 
tive, 

f{q) = ^ and .9(5+) = .g(g_) = 0. (157) 

In the second equality we used Eq. ()155p . which yields 
g{q+) = g{q-), and the fact that g{q) being defined up 
to an additive constant we can choose g{q-) = 0. Next, 
substituting into the first Eq. (|150p . integrating by parts 
and derivating once, we obtain 



q- <q<q+ 



r — q 
t 



1+ 



Aq'CuoM.q')9{q'), (158) 



where we used the first relation ([5^ . Then, we can devise 
a systematic procedure to solve Eq. (|158p when d — n \s 
an odd integer, that is for 



d^n + l + 2i with ^ e N. 



(159) 



Indeed, from Eqs. ([32|) and ([27| we have, with a normal- 
ization factor D, 

1*00 

Cu„AQi,q2)^D dfcfc"$fc((7i)$fc((72), (160) 
Jo 

where we introduced the eigenfunctions on [0,cx)[ of the 
linear operator C, 

q>0: ^fc(g) - /:.$fc-fc2$fe, (161) 



with 



C = - 



[kqy/^ 



d^ 1 - d d d-l 



dq q dq q^ 



(162) 



On the other hand, we note from standard properties of 
Bessel functions (i.e. Hankel transforms) that 



dfcfc''-i$fe(gi)$fe(g2) =gr''<5D(gi -92). (163) 



Therefore, noting C'^ the adjoint of £, 
^1 d^ d - 1 d 



dg^ g dg 



(164) 



we have when condition (I159P is satisfied, for g- < g < 
9+, 



dg'a„„(g,g') ic^' 



jd-ir-q 



Dt 



1+ r°° r — n' 

dg' / dkk^^M g""'— ^ + b.t. 

q- Jo ^ 



r — q 



b.t. 



(165) 
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where we used Eas. (|161l) . (I163p . and "b.t." stands for 
boundary terms at q' ~ q± generated by the integrations 
by parts over q' . Thus, we obtain the solution of Eq. (|158p 
as 

5(q)-/:V-i!:__9+b.t. (166) 

where the boundary terms are of the Dirac type, such as 
SD{q~q±) and its derivatives, locaHzed on the boundaries 
q± . Using the relation (|159p this yields 

^-0: .g(g)(xg"^, £ > 1 g{q) ^ ^ +h.t. 

(167) 

Then, the kernel f{q) can be obtained from Eq. (|157p . 
However, from Ea. p49p we can see that g{q) directly 
gives the velocity profile of the saddle-point as 

r<i+ 

uoq= dq' CuoA<l,<l')9i<l'), (168) 

Jq- 

which follows the linear slope (|158p in the interval 
[q-,q+], while the action (|153p writes 

S^Vrir^"^ j'\qg{q)'-^. (169) 

Since the second equality (|157p has already fulfilled the 
constraint p55p . associated with the minimization with 
respect to the parabola height c, we only need to mini- 
mize S over g_ and g+ to complete the derivation of the 
saddle-point and of the rare-event tails of the distribu- 
tion V{r]r) (after we check that this minimum does not 
give rise to other shocks). In fact, from Eq. (|169p . in the 
present case the analog of Eqs. ([M)) and (|148p reads as 

VM - e-i^'- '^''^^'''>'^. (170) 

Thus, for overdensities q-^- is defined as a function of 77,. 
from the first Eq. (jl56p . and we only need to minimize 
S over g_ to determine the saddle-point and the high- 
density tail (|170p . For underdensities g_ is set by the 
second Eq. p56p and we must minimize S over g+ to ob- 
tain the low-density tail. 

Again, this approach only holds for the rare-event tails, 
in both limits of large scale/early time at fixed density, 
and of extreme density at fixed scale and time, where 
the probability (|170p is much smaller than unity. In such 
regimes, the expression (|170l) gives the asymptotic tail of 
the probability distribution at leading order. 

5. Tails of the velocity divergence distribution 

The method described in the previous section can also 
be applied to the distribution of the velocity divergence, 
7^(8^.). As for the quasi- linear regime studied in sec- 
tion IIVBI the only difference as we go from the overden- 
sity to the velocity divergence is to replace rjr by Qr in the 



action S. In particular, we recover the same saddle-point 
and the same rare-event tail (|170p . and we only need to 
specify the relation between Or and the Lagrangian radii 
q- and g+. Applying the discussion below Eq. (jl56p to 
Qr, which is given by Eq. (jl04p for regular points, we now 
write 

e,. > : Or = - 1) , (171) 

< : Or = - 1) . (172) 

This gives back the relation (|12ip . so that the tails of 
the distribution ^(Gr-) can again be obtained (at leading 
order) from the tails of 'P{rjr) by substituting the second 
Eq.linil). 



6. Case n = G, d = 1: white-noise initial velocity 

Let us describe how this procedure works for the case 
{n = 0, d = 1} (whence ^ = in Eq. p^ V where the 
variance of the linear density contrast is actually 
divergent, so that shocks must always be taken into ac- 
count. Using the initial velocity correlation given in Ta- 
ble |T1 we immediately obtain the solution of Eq. (|158p . 

q-<q<q+: g{q) = g[q±) = 0, (173) 

which gives the hnear velocity profile 

q G [q-,q+] ■ uoq = —j—, q i [q-,q+] ■ u^q = 0. 

(174) 

Note that g{q) is singular (discontinuous) at q±, since 
9{q±) = 0; which gives Dirac terms Soiq — q±) for the 
kernel f{q). Then, the action (jl69p writes 

2 

^-yir+^-^T ~rf + {r^q^ f] . (175) 

For overdensities, rjr > 1, the upper boundary q^ is given 
by p56p . (7+ = rr]r, whereas g_ is determined by mini- 
mizing S. This gives q- = r, since > r is excluded 
as it would give further shocks over the range [r, (;_]: 
we must check that the profile f/jor does not cross the 
parabola Vr.c outside of [q-,q+], that is, that the veloc- 
ity profile does not create a larger shock. For the simple 
profile (|174p we do not need to consider ^or to check 
that no shocks appear beyond [g_, Thus, the system 
is motionless over [0,r[ and -|-cx)[, and particles in 
the range have the linear initial profile (|174p and 

simultaneously merge at radius r at time t. Note that 
there appears a rarefaction interval (empty region) over 
]r, (7+[ as the initial velocity is discontinuous at g+. This 
is due to the large power at high k of the initial white- 
noise energy spectrum (llip . We can see that there are no 
other shocks over disjoint regions that modify the density 
within radius r at time i, so that we have obtained the 
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true minimum over symmetric initial conditions. Then, 
Eq. pTO)) gives 



>1: VM - e-^^^^-^^^^^' \ (176) 



For underdensities, we obtain by a similar reasoning g_ 
rrjr, q+ = r, and 



rjr < I ■ 'P(??r) ^ e" 



(177) 



From Eq. (|12ip we obtain at once for the velocity diver- 
gence Or, which is also the dimensionless velocity incre- 
ment (I22|, the tails 



7'(e.)^e-'-'l«'-l'/(«'). 



(178) 



Since the case {n ~ 0,d = 1} of white-noise one- 
dimensional initial velocity can actually be solved (isl . 
[30| . we can compare the results (|176p - (jl77l) with the ex- 
act distribution Vijir)- Using the notations of [s^l, it 
is known to display the asymptotic behaviors at large 
scales, 



X>1, |r;^-l|>X-\ 7y^>X-3 



and at small scales. 



X < 1, 77^ > X-'^ : r{Tjr) e 



-uJiXrir-X-'rjf./12 



(179) 



, (180) 



where we did not write power-law prefactors, and —uJi is 
the zero of the Airy function Ai(x) closest to the origin 
{uji ~ 2.338). Here X is the dimensionless length of the 
interval [— r, r] of radius r, whence of size x = 2r, 



X = 



2r 



2r 



(2£li2)l/3 (4i2)l/3^ 



v-3 

hence = (181) 

12 6^2' ^ ' 



since the normalization used in the present paper cor- 
responds to D = 2 [6l|. Thus, we can check that for 
large overdensities our saddle-point result (|176p agrees 
with the exact results (|179p - (|180l) at leading order, at 
both large and small scales. Of course, this only applies 
to the rare-event tail, which is repelled to larger densi- 
ties, rjr ^ l/A', at small scales in the highly nonlinear 
regime. For large underdensities, we also recover the ex- 
act result p79p at leading order, that applies to large 
scales. It cannot give the low-density tail in the highly 
nonlinear regime because this no longer corresponds to 
rare events. Indeed, as described in [30], at low densities 
the distribution V{rir) shows an inverse square root tail, 
oc 1/y^, and a Dirac contribution, Soii^r), that both 
have a weight, of order q-<^i^-^^ /^^ at large scales, that 
becomes of order unity at small scales. In fact, on small 
scales most cells of radius r are actually empty, so that 
there no longer exists a rare-underdensities tail. Note 
that this feature can actually be seen on the saddle-point 
result (|177p . as we can see that for r < t^/^ the exponen- 
tial becomes of order unity for 77^ = 0. This shows that 



empty or almost empty regions are no longer rare, and 
that the distribution Virir) over this range cannot be ob- 
tained by a saddle-point approach of the type described 
in this article. 

Finally, it is interesting to note that Eq. (|176p agrees 
with the behavior that would be obtained by a naive 
extension of Eq. (|133p to {n = 0, d = 1}, even though 
the derivation of Eq. (|132p does not apply to this case 
(the variance cr|^^ is even divergent). On the other 
hand, for underdensities Eq. (|134p would give ln7-'(?7^) ^ 
—t^'^r^rj,. for 77,. 0. This shows at once that 
this cannot give the low-density part of the probabil- 
ity distribution, since we do not find a rare-event tail 
{hiV{rir) —00) but a probability of order unity, which 
a saddle-point approach cannot describe. As discussed 
above, this is not a failure of Eq. p34p . since there is no 
rare low-density tail as empty regions occur with a finite 
probability, that goes to unity at small scales. 

Of course, the discussion above also applies to the dis- 
tribution V{Qr)- Thus, the tail (|178p agrees with the 
exact result for large positive Or in all regimes, and for 
negative Qr in the quasi-linear regime, while on small 
scales, in the highly nonlinear regime, cells with O^ ~ — 1 
(associated with almost empty domains) are no longer 
rare and cannot be described by the method used here. 
Again, the scalings obtained in the exponential (|178p 
agree with a naive extension of Eq. (|135p while the ex- 
tension of Eq. (|137p correctly signals the absence of rare 
low- Or tail. 



7. Case n = 0, d = 3 

We now consider the case {n = Q,d — 3}, which gives 
e = 1 in Eq. p^ . From Eq. p^ the regular part of 
g(q) is proportional to r/t, and we find for the solution of 
Eq. pSSp with the normalization of Cug^ given in TableHl 

9W ^ + ^D{q - q-) + Soiq ~ q+). 

(182) 

Here the Dirac terms should be understood as Soiq — 
{q± =F e)] with e ^ O"*" (i.e. they have an integral weight 
of unity within [g_,g+]). This gives the linear velocity 
profile 

q{r - 9-) ^ r 1 r-q 
q<q-: Uoq = , q £ [q-,q+\ ■ "Og = 



q> q+: Uog 
and the action 



tq- 
ql(r-q+) 
tq^ 



S = yvr + ^ (4r2g+ - 6rql + 3ql - r\^) 



t ' 
(183) 

(184) 



For overdensities, 77^ > 1, the minimization over q_ 
gives q_ = q+- Indeed, contrary to the previous case, 
{n = 0,d — 1}, it is now possible to have q^ > r without 
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building a larger shock, as already seen from Fig. [2l since 
we actually recover the saddle-point of section IV Al with- 
out shocks, as — (7+ (i.e. an isolated contact point 
between the parabola Vr,c{q) and ipoq)- Then, Ea. (|170p 
writes as V{rir) ^ e~^''^^^ '^-^l ^)^/(2t^)^ j-^ agreement 
with Eq. fTUD)) and Table El and we recover the tail (fT^ 
and Table rvTITl Indeed, the constraint in (|132p is satis- 
fied, so that we already knew that we had to recover that 
regular saddle-point. 

For underdensities, 77^ < 1, the minimization over q+ 
gives q+ = 2r/3, so that we only have a shock (i.e. 
g_ < q+) for g_ < 2r/3, that is for low densities be- 
low = (2/3)'^. This agrees with the discussion in sec- 
tion [IVXJ] and Table [TTIl where we found that the regu- 
lar saddle-point ([7S)l only develops a shock after a finite 
time, that is below a nonzero density contrast thresh- 
old. In the quasi- linear limit, r^/t'^ ^ 1, where the 
range ry* < rjr < 1 already corresponds to large den- 
sity fluctuations, we can also use the method described 
in section IIVI and we obtain the result of Tabic |Vl The 
analysis described above from the action ()184|) provides 
the density threshold r|^, = (2/3)'^ written in that Table. 
From the relation p2ip , this also gives the velocity diver- 
gence threshold O* — —1, above which the quasi- linear 
distribution 7^(0^) is given by Table rvTIl 

For larger underdensities, < rir < ry*, we have q- < 
q+ and we must use the action (|184p that takes shocks 
into account, since we can check that the profile (|183p is 
valid (there are no other shocks that modify the density 
within radius r). Then, Eq. (|170p gives 



0<7^<(2/3)^: InVM 



6^2 V 9 



(185) 



quasi- linear and highly nonlinear regimes. For low Or, 
Eq. (fT85l) gives 



-3<er<-l: InVier 



1 



18t2 13 



- + e^ . (186) 



Again, this provides the very low-Qr tail in the quasi- 
linear regime, which disappears in the highly nonlinear 
regime where there is no longer a rare-event \ow-Qr tail, 
and this behavior can also be seen in the naive extension 
of Eg.lfnTl). 



Case n - 



-1, d ■- 



We now turn to the case {n = — 1, d = 2}. As seen from 
Table IVlIIl shocks should only appear for underdensities, 
as in the previous case {n = 0,c? = 3}. The extension of 
Eq. ()134p gives a vanishing power of rjr , so we can expect 
a logarithmic dependence on rjr (or a finite asymptote) 
for In 7^(77,.) at low densities. We again have ^ = 1 in 
Eq. (|159p . so that the regular part of g{q) is obtained 
from Eq. (|167p as cx r/{tq), and we find 

^('^^ ^ i^+^'^^^'^"'^-)+^^^r^^^('^"'^+)' ^^^^^ 

and 

q<q-: uoq = , qe[q-,q+\ - uoq = 



q> q+: U09 



tq- 

q+{r-q+) 

tq 



t • 
(188) 



while the action writes 



Of course, we can check that at point r\r = ?7* Ea. (|185p 
is equal to the result (|100p . shown in Table fVl which is 
provided by the regular saddle-point. We can see that 
below this threshold the dependence on 77^ of V{rir) is 
modified by shocks. Thus, Eq. (|185p provides the very low 
density tail of 'P{rjr) in the quasi-linear limit, /t^ ^ 1. 

In the nonlinear regime, r^/t^ < 1, the result p85)) 
becomes of order unity over the range Q < rjr < ??* . Note 
that this agrees with the naive extension of Eq. (|134p . 
which also yields the correct exponents of r, t and rjr. 
Then, as for the case {n = 0, = 1} studied in the pre- 
vious section IV B 6) there is no rare underdensities tail, 
and empty or almost empty regions are not rare. More 
precisely, there is no exponential decay of V{rjr) at low 
rjr, but power-law prefactors associated with subleading 
order terms may give either a power-law growth or falloff 
at rjr — !■ —d. However, the precise behavior of the dis- 
tribution V{rjr) for rjr ^ cannot be derived through a 
saddle-point method since there is no rare tail and one 
should take into account many possible initial configura- 
tions. 

The previous results directly extend to the distribution 
V{&r)- Thus, for large 8^ we recover the saddle-point 
of section IV Al and the tail (|135p . which applies to both 



^ = + y ^ + ~ ^'''^+ + ^'^+) ■ ^^^^^ 

As expected, for overdensities we recover g_ = q^ (i.e. 
the regular saddle-point without shock) and 'P{rjr) ~ 
e-'-'(v^-i)V(2t")^ in agreement with Eq.fTUD)) an d Ta- 
bleEl and we also recover the tail (fT32)) and Table I Villi 

For underdensities we obtain q^ = r/2, so that we 
only have a shock below rj^ = 1/4, which provides the 
density threshold written in Table |Vl Thus, as for the 
case {n = 0, (i = 3} and in agreement with section llV A 21 
in the quasi-linear regime for rare underdensities in the 
range 77* < 77^ < 1 the distribution V{rjr) is obtained 
from the method described in section HVl which gives the 
result of Table IVl In terms of the velocity divergence Qr, 
this also provides the threshold 8* = — 1 of Table IVTIl 
and above this threshold the distribution 7^(8^) is given 
by Eq. (|119p and Table IVlTl in the quasi- linear limit. 

For lower densities, in the range < ry^ < ?7*, we obtain 
from Eq. (fT89l) 



< 7/^ < 1/4 : \nV{fjr) 

whence 'P{rjr) ^ (4?7, 



8t^ 



[l-ln(4r;,)], (190) 



(191) 
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n 


d 


In 75(77,.) for rjr > 1 


In 75(77,.) for 77r < 1 


In 75(9,.) for e,. > 


In 75(9,.) for 9,. < 





1 













3 


no shock 


{l-nl") for 7;. < A 


no shock 


^(i + 9.) for9.<-l 


-1 


2 


no shock 


[l-ln(47,0] forr,,. < \ 


no shock 


-^[l-21n(2 + 90] for9r<-l 



TABLE IX: The tails of the distributions V{rjr) and 75(9r) for the initial conditions of Table H] in cases where shocks appear. 
These results also hold for the quasi-linear limit, t ^ or r ^ 00, at fixed rjr or 9r, but only below the thresholds given in this 
Table for the last two rows (higher densities and velocity divergences are described by Tables IVl I VIII) . For {n = 0, d = 1}, these 
results hold in the quasi-linear limit for any r/r 7^ 1 and 9r 7^ 0. In the highly nonlinear limit, r ^ or t — > 00, the rare-event 
tails at low 77^ and 9r disappear as low densities and velocity divergences are no longer rare (ln7' in this Table becomes of 
order unity and these formulae are no longer valid). Cases marked as "no shock" correspond to saddle-points that do not give 
rise to shocks, so that the results of Tables FVl I VIII and I Villi are valid. 



Again, at the transition 77^ — rj^, Eq. (|190p is equal to 
Eg. poop shown in Table IVl Thus, we obtain as expected 
a logarithmic dependence over 77,. for ln7-'(?7r), in agree- 
ment with (|134p . In the quasi-linear regime, r^/i^ ^ 1, 
Eq. (|19ip means that the low-density tail has a power- 
law behavior V{rir) ~ 77", with an exponent a ^ j (8i^) 
that grows at large scales and early times, so that the 
low-density fallofF is increasingly sharp. However, be- 
cause there could be a power-law prefactor in Ea. p9ip 
due to sub-leading corrections to the steepest-descent ap- 
proximations, this is unlikely to give the exact exponent 
a but only its behavior at large r and small t. Then, in 
the nonlinear regime, jt^ <^ 1, Eq. (|19ip is not sufficient 
to give the behavior of P(7?r) for rjr 0, as these pref- 
actors may either give a positive or negative exponent. 
This limiting configuration between the cases n > — 1, 
where empty or almost empty regions have a finite prob- 
ability at small scales, and n < —1, where low densities 
exhibit an exponential tail of the form (|134p . requires a 
finer analysis in the nonlinear regime. 

For large velocity divergence Qr we recover the tail 
(I135P associated with the regular saddle-point while for 
low e,. Eqs.dHOl) and p2T|) yield 

- 2<er<-l : In 7^(9,.) - — [1-2 ln(2 + er)], (192) 

whence ^'(9.) ~ (2 + 9.)'''/'^*'^ e-'^'/(«*').(193) 

This only gives the behavior at low 9r in the quasi-linear 
regime, r'^/P ^ 1, as in the nonlinear regime power- 
law prefactors may lead either to a growth or decay of 
75 (9r), but in both cases there is no rare-event tail (i.e. 
no exponential falloff). 

9. Summary for low integer n and d 

We summarize in Table IIXI the results obtained from 
the approach developed in the previous sections for the 
tails of the distributions V{rir) and 7^(9^), for the initial 



n 


d 


ln[n(m)] for 7n ^ 00 





1 


-7nV(48tVo) 





3 


-3m/(87rtVo) 


-1 


2 


-m/(27rtVo) 


-2 


1 


-77i/(4tVo) 


-2 


3 





TABLE X: Large-mass tail of the mass function n(m) of 
point-like objects, for the initial conditions of Table U] from 
Ea. (IT96|l . 



conditions of Table [T] where shocks cannot be neglected. 
This complements the Table IVIIII that applies to cases 
where the saddle-point does not form shocks. 

These rare-event tails apply to the large fluctuation 
limits, r]r 00 and rj^ — > 0, or 9^ — > 00 and 9^ —d, 
at fixed time and scale, whatever the value of the variance 
ai or afi, . 

They also apply to the quasi-linear limit, 175^^ ^ or 
fi/'or ~* 0; that is at early times or large scales, at fixed 
r]r and 9^, below some finite thresholds rj^ and 9* in the 
two cases {n = 0, d = 3} and {n — —l,d = 2}, and for 
any rjr ^ 1 and 9r 7^ in the case {n ^ 0, d = 1}. 

For these three cases, in the highly nonlinear regime, 
r —^ or t 00, the rare-event tail at low rjr and 9^. 
disappears as these results predict that InT' becomes of 
order unity. Then, low densities and velocity divergences 
are no longer rare (but the probability distribution might 
still decay as a power law) and cannot be described by a 
saddle-point approach. 
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C. Mass function of point-like singularities 



D. Pre-shock singularities 



As the density and velocity fields evolve through the 
nonlinear Burgers dynamics, starting from the scale-free 
initial conditions the system displays an intricate 
self-similar progression from smaller to larger scales. In 
particular, collisions between particles create discontinu- 
ities (shocks) of dimension d—1, d—2, .., down to 0, lower 
dimensional objects arising at the intersection of higher- 
dimension structures. For instance, if d = 3, once parti- 
cles have formed a two-dimensional sheet of finite surface 
density, orthogonal to the direction of the largest eigen- 
value of the initial tidal tensor, they keep moving within 
this surface and form critical lines and nodes. Then, the 
typical distance between such objects increases as L(t), 
as in ([H]), and their mass grows accordingly. The mass 
and the overdensity within radius r about a point x, con- 
tained in such a £)-dimensional structure, scale as 







m{< r) ^ fJ-r^, 
to(< r) 



Vr 



„D-d 



(194) 
(195) 



Thus, at small scales we can see that very large densities 
are associated with the lowest-dimension objects, D = 0, 
and the contribution of these point-like masses to the 
probability distribution V{rir) reads as 

r — !■ 0, rjr oo : V{rir)drir ^ Vn{'m)dm, (196) 

where V is the volume of the sphere of radius r and 
nlm) is the mass function of point-like masses, that 
is, n(m)dmdx is the mean number of such objects of 
0— dimension, with a mass in the range [m,m + dm], 
within the volume element dx. Then, from Eq. ()133|) we 
obtain for the high-mass tail 



ln[n(rn)] oc — - 



(197) 



Indeed, we have seen in scction lVBI that the scaling (|133p . 
that was derived in section IV A II for n < d — 3, actually 
extends to the full range —3 < n < 1, but the propor- 
tionality factor is no longer set by Eq. (|132p . However, 
in the range n < c? — 3, this numerical factor is given by 
Eq. ()132p , while for d — 3 < n < 1 it can be obtained from 
the analysis described in section IVBl and from Table HXl 
for the associated integer values n = and d—1. We 
show our results in Table |X] for the high-mass tail of the 
mass function n(m) of point-like objects, for the initial 
conditions of Table HI 

Of course, as for the density and velocity distributions, 
these results agree with the exact expressions that can 
be obtained in the two cases jji = 0,d — 1} and {n — 
-2, = 1} [ii, [H, m 113, H S i^. For more general 
cases, the scaling (|197p was already conjectured in [ll], 
[T^ . from numerical simulations and heuristic arguments, 
and proved in 19] for — 1 < n < 1 with d = I (with 
upper and lower bounds for the proportionality factor). 



Before we conclude, we should add a few comments 
on the comparison of this work with studies of pre- 
shock singularities As shown in [s^], for smooth 
initial conditions large densities are localized near "kur- 
toparabolic" singularities residing on space-time mani- 
folds of codimension two. They lead to universal density 
tails V{r]) ~ ry"''''^ in any dimension. In one dimen- 
sion, this corresponds to pre-shocks [H,!!^, that is, when 
a shock forms the Lagrangian potential changes from a 
single extremum to three extrema and at the transition, 
where its second derivative vanishes, one can see through 
a Taylor expansion that the Eulerian density field be- 
haves as x~'^l'^ close to the singularity. Then, the con- 
tribution from the neighborhood of such events (both in 
space and time) yields a power-law tail V[r\) ~ r?"''^^. 
This can also be extended to higher dimensions [1, 113] ■ 
These results apply to the unsmoothed density field for 
smooth initial conditions. By contrast, in the present ar- 
ticle we study the smoothed density and velocity fields, 
that is we always consider the mean density and veloc- 
ity increment over a finite radius r, for non-smooth ini- 
tial conditions described by the power-law power spec- 
tra (|18p . Thus, these are two very different regimes. In 
particular, this explains why we obtain probability dis- 
tributions that depend on both the dimension d and the 
slope n of the initial power spectrum (over the range 
—3 < 71 < 1), rather than universal tails. We can note 
from Eq. (|133p that in the regime studied here the proba- 
bility distribution Vijir) shows a large-density exponen- 
tial tail with a characteristic cutoff T]r ^ 'r~'^ that goes 
to infinity as r — > 0. Then, at very small scales (i.e. 
in the highly nonlinear regime) an intermediate power- 
law regime can develop below this upper cutoff. How- 
ever, this power law is not universal either, since the 
exact results obtained in [H, [s^] show that for d = 1 we 
have 'P{rir) ^ vir^^'^ if n = —2, see also Ea. (|103p . and 
V{rir) ^ Vr a n = 0. 



We can note that for the forced Burgers equation 
similar universal power-law tails can be obtained us- 
ing instanton methods (i.e. looking for relevant saddle- 
points, that correspond to shocks) [54], although there is 
some debate on the exact value of the exponent, which 
might depend on the influence of the boundary condi- 
tions [5^ [13] ■ Again, these results consider a smooth 
forcing so that the exponent is set by the dynamics of a 
single shock and is universal. For singular forcing (i.e. 
with significant power at high wavenumbers) one might 
obtain non-universal results for the density and velocity 
increments over finite radius r, in a fashion similar to 
the free case studied here. However, this goes beyond 
the scope of this article. 
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VI. CONCLUSION 

We have studied in this article some asymptotic prop- 
erties of decaying Burgers turbulence in d dimensions. 
Focussing on the case of random Gaussian initial veloc- 
ities and a uniform initial density, we considered power- 
law initial energy spectra such that the evolution is self- 
similar. Thus, the system displays a hierarchical evo- 
lution and the integral scale of turbulence, L{t), that 
is generated by the Burgers dynamics and separates 
the large-scale quasi-linear regime from the small-scale 
highly nonlinear regime, grows with time as a power law, 
L(t) oc Then, in order to take advantage of the 

statistical homogeneity and isotropy of the system (once 
we have taken care of the infrared divergence if n < —1), 
we have defined the spherical quantities, rjr and 6^, that 
are the overdensity and the velocity increment over a 
sphere of radius r. This allows to preserve the statistical 
isotropy of the problem and to handle the case of large 
dimensions d> 1. 

We have first recalled how such a nonlinear dynam- 
ics can be studied through standard perturbative expan- 
sions. Here this corresponds to expansions over powers 
of time, or equivalently over powers of the initial velocity 
fluctuations. This approach is quite flexible, as it does 
not require any symmetry, but it becomes very heavy 
at high orders. It can be somewhat simplified when one 
focusses on spherically symmetric quantities such as r/r 
and Or, defined through a real-space top-hat filter, but it 
remains cumbersome for arbitrary dimensions. We have 
pointed out that from a perturbative point of view the 
Burgers dynamics in the inviscid limit is equivalent to 
the Zeldovich dynamics. This means that shocks are not 
taken into account and require non-perturbative meth- 
ods. 

Next, we have described how to derive the asymptotic 
probability distributions V{r]r) and reached in the 

quasi-linear regime from a saddle-point approximation. 
This method allows to obtain at once the asymptotic 
cumulant generating function (p{y), the Taylor expan- 
sion of which provides the leading-order term for each 
cumulant {r]^)c, that would be obtained from the previ- 
ous perturbative expansion truncated at order p — 1. In 
addition, the generating function ip{y) being obtained di- 
rectly trough a steepest-descent computation, one can go 
beyond its apparent singularities (associated with large 
high-order cumulants and slowly decaying tails for the 
probability distributions) and make sense of possible sec- 
ondary branches, that are found when this function ap- 
pears to be multivalued. This approach takes advan- 
tage of the spherical symmetry of the observables rjr and 
Or to reduce the problem to a one-dimensional system, 
as the relevant saddle-point (instanton) is also spheri- 
cally symmetric. This allows to derive simple results for 
arbitrary dimension d, provided this instanton has not 
formed shocks yet. 

Then, from the radial profile of this saddle-point, we 
have found that these results only apply to the range of 



initial energy spectrum index —3 < n < d — 3 (within 
— 3 < n < 1), and only above a nonzero imderdensity if 
—2 < n < d — 3. For n > d — 2 the quasi- linear regime 
does not really exist. More precisely, the overdensity and 
velocity divergence rjr and 9^. are already divergent at lin- 
ear order and their distributions do not converge towards 
a Gaussian at early time or large scale, in spite of the 
Gaussianity of the initial conditions. Thus, the system is 
always dominated by shocks. For d — 3<n<d — 2 the 
linear theory is well defined, so that one recovers Gaus- 
sian distributions at very early times or large scales, but 
the saddle-point forms shocks as soon ast > 0. Then, the 
qualitative results obtained from this steepest-descent 
approach should remain valid and still provide a reason- 
able quantitative approximation, as shocks appear over a 
limited range of radii, but they are expected to be mod- 
ified by prefactors of order unity. 

Thus, in order to describe the cases d— 3<n< 1, as 
well as very large underdensities for — 2 < n < d — 3, it 
is necessary to take into account shocks. We have shown 
how to modify this saddle-point method, taking advan- 
tage of the geometrical interpretation of the Hopf-Cole 
solution in terms of first-contact paraboloids, to handle 
these cases. This allows us to find out the instantons, 
which contain shocks, that provide the leading-order be- 
havior of the rare-event tails of T'(?7r) and 7'(0,.). Fo- 
cussing on some low integer values of n and d, where sim- 
ple explicit results can be derived, we have obtained the 
asymptotic tails of these probability distributions, at any 
finite time and scale, for the cases {n, d} = {0, 1}, {0, 3} 
and {—1,2}. We note that the scahngs actually agree 
with a naive extension of those obtained from the regu- 
lar saddle-point computation. This also gives the high- 
mass tail of the mass function of point-like singularities 
(i.e. Dirac peaks in the density field, which correspond 
to shock strengths in the one-dimensional case). 

Then, we find that for n < —1 the very low density 
tail shows an exponential cutoff of the form e"''^"^^''^, 
whereas for n > — 1 there is no exponential falloff (but 
there could be a power-law decline). For the latter cases, 
in the quasi-linear regime, this part of the probability dis- 
tribution V{rir) corresponds to extremely rare underden- 
sities and has a negligible weight, and around moderate 
fluctuations, jry^ — 1| <C 1, the distribution shows a falloff 
on both sides of the mean (77^) = 1. In the highly nonlin- 
ear regime, this intermediate low-density regime disap- 
pears and low-density (and empty) regions are no longer 
rare. Then, one needs another method to describe the 
low-density part of the probability distribution at small 
scales. 

Throughout this article, we have checked that our re- 
sults agree with the two one-dimensional cases of white- 
noise initial velocity (n = 0) and Brownian initial velocity 
(n = —2), where many exact results are known, thanks 
to the Markovian properties shared by both cases, which 
allow a derivation of explicit formulae through specific 
methods. Note that these two cases are representative 
of the two classes of initial conditions, — 1 < n < 1 and 
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—3 < n < — 1, where the initial velocity is dominated by 
small/long wavelengths and which exhibit the different 
behaviors discussed above. Hence they provide a good 
check of the general methods presented in this article. In 
addition to the interest of the asymptotic behaviors ob- 
tained here, we can hope that they could serve as a bench- 



mark to test other approximation schemes, devised to 
study additional quantities such as typical events. More- 
over, since the approach developed in this paper is rather 
general - for instance it was already applied to the grav- 
itational dynamics - it may also prove useful for other 
systems. 
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